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Chapter 1 

Symmetries in Field Theory. 



1.1 Introducing Q.F.T. and Symmetries. 

The generating functional for a quantum field theory containing a real scalar field is given by the following 

Z = j D(l)expi j dPxC{(l)) (1.1) 



where 



The correlation functions in this scalar field theory are given by 



C {ct>) = ^V'^'d^cl^d,^ - -m^4>^ -V{cl>). (1.2) 



G{xi,...,Xn) = ^ j D(t)(t){xi) ...<t){xn) e^vi jd^x£{(j)). (1.3) 

Suppose the potential V (0) has the form V {(j)) = gcj)^. In this instance the Lagrangian has a synmietry cj) — > —cf) as 
L {(j)) is an even function of ^. Before continuing consider what is meant by the measure of the integration, J D(p. 

D<I> = ]1 d<^(x) (1.4) 

^ J —CO 

One should consider x as being a label for a variable, the variable being the value of the field at the space-time point x . 
Hence the product of the integrals of the values of the fields at every point in space-time is the integral over all possible 
field configurations. So what happens to the integral over all possible field configurations under ^ — > ^' = —cj) ? 

jDc^ = n.rr#(^) 

= n./+"r-#'(^) ,151 
= 

Nothing, the measure is invariant under (p ^ (j)' = —(p. Now as the measure of the integration and the action are both 
invariant under this transformation of the fields the whole generating functional is therefore invariant under 0^0' = 
-(j), Z ^ Z. 

The effects of this transformation on the correlation functions are not so trivial though, 



G{xi,...,Xn) = ^ J D(j)(l){xi)...(l){xn) expijd^ X £{(!)) 

= ...(/.' (x„) expi/d^:r £(</.') 



(1-6) 
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a factor of (—1)" is acquired from the string of fields before the exponent (j) {x\) ...(j) (a;„) = (—1)" (j)' {xi) ...(f)' (a;„). 
From a purely mathematical perspective one sees that, 

-| j D(^(j>{xi)...<^{xn) expi j d^xC{(/>) = ^ j 0' (xi) ...0' (x„) expi j d^x£(0') (1.7) 

as the field is just a dummy integration variable which we are integrating everywhere from +oo to — oo and because 
Z ^ Z under the transformation. Hence, 

G(xi,...,X2) = /-D(?!)(/)(xi) ...(/) (x„) expi/d^x £((/)) 

= (-l)"|/i?0'.^'(xi)...<^'(x„) expi/d^x£(.^') 
= (-iri/£>00(xi)...0(x„) expi/d^x£(.^) 

= (-1)"G(X1,...,X2) 

For this statement to be self consistent n must be even, i.e. no odd Green's functions. This result was seen in our 
perturbative treatment of QFT in other courses where it was not so much proved as shown. The point to take away here 
is that one can make important deductions in field theory through symmetry considerations and without resorting to 
perturbation theory. Also, we have derived this result in an arbitrary number of space-time dimensions D. 

We can still however raise questions about the vaHdity of this result. One obvious point is that this result is restricted 
to potentials which are even functions of (}>. Also, do quantum loop corrections affect our result? Finally note that the 
symmetry transformation law used to derive our result was discrete. 

Consider now continuous symmetries. Let us take a similar Lagrangian to the one above but this time with two real 
scalar fields (pi and (f)2. 

1 1 

C{(Puh) = ^t''d^^id,(t>2 + ^t''d^(t>id,(t>2 - -Y {<t>l + <t>D -V{<t>iA2) . (1.9) 

The integration over all field configurations is obviously / Dcpi J D(p2. We will now employ a cormnon streamUned 
notation for these situations. We re-express our fields as follows 

= 75(<^i-#2) \(P^ J V2\l -i J {(^2 J 

Note that the dagger here is not that of canonical quantization, 0^ is not an annihilation operator! If that was the 
case we would only have one type of creation operator and one type of annihilation operator in decomposing </> and 
this is not the case! has its own creation operator and its own annihilation operator and ^2 also has its own 
creation operator and its own annihilation operator i.e. there would be two distinct creation operators and two distinct 
annihilation operators. In the new compact notation we have, 

mass term = — m^0^^ (1.11) 

kinetic term = 1,,^. a. (^) } + i,- (^) z^. (^) } 

= ^'n^"' { \ (2a^(/)ta^0 -h 2d^(l)d^(j)^) } ^ Use r]'^'' fj, ^ u symmetry (1.12) 
= v'"'d^,(P^d^(P 

The action is then, 

5[(/),(/)t] = J d^x {ri''''d^,(P^d,(j,-m^(t)U-V{^,(t>^)}. (1.13) 
Under a global transformation e'" e i7 (1) of the fields: ^ — > 0' = e'"(^. The kinetic and mass terms in the Lagrangian 



1.1. Introducing Q.F.T. and Symmetries. 



3 



are invariant: 

n''''d^4>^d,ct> = ri^^^d^ (0'te+-) (e-*"<^') = r/'^^ {d^cj,'^ d,<P') e+^^e"*" = tj'^" {d^4>'^ d,4>') (1.14) 

m^t,^ = mV^e+'«e-*«(^' = w? cj)"^ cj)' . (1.15) 



Hence to make the whole Lagrangian invariant we need to have a potential V [(j), (j)') which is invariant. This is clearly 
the case if the potential is a function products of fields of the form t^^t^. Hence we have that the action, 

S[(j),(j)^ = J d^x {r]^"'d^(j)''d,(l)-m^(f>U-V{(l)U)} (1-16) 

possesses a global U {1} invariance. 

What about the measure of the integration? 

D(f>i Dh=l[ d(/)i (x) / d(l>2 (x) . (1.17) 

J ™ J —oo J — oo 



i)' = e*^0 = {(j)\cosa — (j)2sina) + i {(l)\sina + (f)2Cosa) = (j)'^ + i4>2 



^/f _ g «a^t = [(j)icosa — (j}2sina) — i {(j)isina + (j)2Cosa) = (j)'^ — icf''- 



(1.18) 



\ / cosa -sina \ ( (t>i \ / IItIjtW'/'iI t/'^i 



sina cosa I \ <!>■ 



J\ (1.19) 



/!■ r-+oo r-+oo / r-+oo /.+00 \ 

Dcp[ d0'i (x) / d0^ (x) = [] / / (Det J (x)) d0i (x) d02 (x) (1.20) 

./ 2: J —00 J — 00 ^ \J — 00 J — CO J 

Note that the Jacobian in the integral is equal to one for all x (Det J (a;) = 1)! 

//• r+oo r+00 /■+00 /•+00 r r 

D<l>[ / £>02 = n / '^'^i / '^'^s (^) = n / (x) / d<^2 (a;) = / i'-Ai / £>02 

3. J —00 00 ^ J —00 00 

(1.21) 

i.e. the measure of the path integral is also unchanged by the global U (1) transformation of the fields. What happens 
to the correlation functions? 

G{xi,...,Xn,Xn+l,...,Xn+m) = ^ j Df, {(j), (j)^) (j) (xi) ...(/) {Xn) 'l)Hxn+l) ■■■<l>Hxn+m) ^W^S [<l>, <l>^ (1.22) 

In the above Green's function the variables xi, Xn apply to the (f)'s and Xn+i, ■.■,Xn+m apply to the ^^'s (note that 
/ Dfj^ denotes the integration over the all possible configurations of the fields (j) and (j)'). Now consider the 

effect of the U (1) transformation (f> (x) 4>' (x) = e+*"(/) (x) and 0^ (x) (f>'^ (x) — e~^°'(j)^ {x). The integration 
measure and the action are unchanged as we reasoned above, the string of fields before the exponent, acquires a phase 
factor e'(™~")"(/)' (xi) ...(j>' {xn) (t>'^ {xn+i) ■■■(t>'^ (a;„+m). We rewrite the Green's function in terms of (j)': 

G{xu...,Xn,Xn+l,...,Xn+m) = e'(™-")"^ / (e"- 0', e+*"(/)'t) Cj)' {Xr) ..4' {Xn)cj>'Hxn+l) ..4'Hxn+m) 

(1.23) 

At the risk of being verbose the measure is: 

J {cj>, <^t) =Jd^ (e— <^', e+-<^'t) = J (^/, ^/f) . (1.24) 
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G{xu-,Xn,Xn+l,...,Xn+^) = c'^'""")"^ / 1?^ (0', </.'t) 0' (xi) (/.'t (a;„+i) ...(/.'t (a;„+„) 

(1.25) 

Note, as before when considering the symmetry (p ^ —(p, one can step back forget any Physics or symmetries going on 
and note that mathematically / (0', (j)'^) (p' (xi) ...(/>' (a;„) (j)''^ {Xn+i) ■■■(p'^ {xn+m) x exp iS [</>', 0''^] is exactly 
the same as ;| / 1?^ (</>, (xi) (a;„) 0^ ...</)^ (x„+m) x exp iS \<p, <p'^\ that is to say both are equal 

to G{xi, ...,Xn,Xn+i, ...,Xn+m)- This IS bccausc the (p and (p^ fields are essentially just "dunrniy" variables for 
the integration and they are all integrated from — oo to +oo everywhere {i.e. the limits of the integration are also 
unchanged - this was implicit in our proof that the measure was invariant and is crucial for the last sentence is to be 
correct). Analogous to the ^ — * —cp case one here has, 

For these two statements to be consistent we must have m = n. Therefore only Green's functions that contain equal 
numbers of the different types of fields (p and (p^ are non vanishing! 
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1.2 Local Symmetries. 

Consider a local U (1) transformation: ^ </>' = e*"*^^-*;/), 0^ ^ 0'^ = e~*"*^^-'0^. Under this transformation the 
measure of the path integral is unchanged as it was in the global case, the only difference is that the transformation of 
variables is different at every point in space-time but is nevertheless always of the form: 

j^^^^f cosaix) -sinaix)\ ^^ ^7) 

\ sin a {x) cos a (x) I 

and hence the Jacobian (Det J (a;)) at every space-time point is always one as before! 

= n.te/-r(DetJW)d0U^)d0'2W) (1-28) 
= jDcl>[jDct>', = jD,{ct>',<f>'^) 



As before if we construct our potential such that it is a function of (p^cp only, then it too must be unchanged by our 
transformation of the fields. 



The mass term is also of the form and hence it too is invariant under these U (1) transformations of the fields. 

Once again our theory is starting to look as though it possesses U {1) invariance, though this time the invariance is 
local. The kinetic term in the action is however definitely not invariant. Whereas before our fields acquired a constant 
phase on application of the ?7 (1) transformation now they have acquired a phase with a space-time dependence hence 
we cannot merely pull the phases out in front of the derivatives in the kinetic term and have them cancel. Instead we 
get the following mess for the kinetic part of the action, 

= fd^x r|^''' (e-*"(^) (9^0' - i {d^a {x)) 0') x e+*"(^) + i {d^a {x)) 0't)) 

= / d^x LKinetic {<!>', 0'^) + / d^x (* (5^<^') d^a {x) - i (a,0't) cp'd^a (x)) 
+ J d^x T?^-^ {d^a (x)) {d,a {x)) cl)'^'^ 

Exploit r]^'^ symmetry in 2nd term (linear in a (x)). 
= J d^ar LKinetic (0', 0'^) + / d^a; {i (9^0') 0'ta,a (a;) - z (a^0't) cj^'d^a (x)) 
+ J d^a; <^ {di,a (x)) {d,a (x)) 0'0'^ 

= Smneuc [0', 0'^] + / d^x ir?^" ((9^0') 0't - (a^0't) 0') (d^a (x)) 
+ J d^x r,^'^' {d^a (a;)) (5,a (x)) 0'0't 

(1.30) 

Now we have to play some games with the term linear in a [x). Integrating by parts we can write it as, 

/d^x in^- ((9^0') 0't - (a^0't) 0') {d,a (x)) = /d^x zr^A'-a, (a (x) ((9^0') 0't - (a^0't) 0')) 

- /d^x 11^^- a (x) ((9^0') 0't - (a^0't) 0') 

The first term J d^x irj^^^diy [a (x) {{d^4>') 0't — (a^0't) 0')) is a Z) — divergence which we can write as a surface 
integral using Gauss's law. The divergence is integrated over all of space-time, so if we used Gauss's law the surface 
over which we would integrate is at infinity. Assuming that the fields are vanishing at infinity limaj^oo {x) = the 
integrand will be zero at all points on the aforementioned surface, hence the integral is zero and we have, 

j d^x iri^^ ((9^0') 0't - (a^0't) 0') {d^a (x)) = -Jd''x irf^a (x) ((9^0') 0't - (a^0't) 0') . (1.32) 
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We shall denote i {dfj,(j)') c^'^ — i {dij,(j)''^) (j)' by j^, so the kinetic part of the action becomes: 

SKmeUc[(t>,(t>^] = S Kinetic [(!>' A'^] -j' ^'^ X a {x) d^j.j^' 

If we take the transformation to be infinitesimal i.e. a {x) very small we can neglect the term quadratic in a {x). So 
overall we have that everything is invariant except for the kinetic part of the action. 

/£>^(</.,<^t) = /D^(<^',</)'t) 



/ . N / . N (1-34) 

SKinetic [<l>, (1)^ = SKinetic [</>', - J d^a; a {x) d^j^ 



Now consider the Green's functions: 



G{xi,...,Xn,Xn+l,-;Xn+m) = j Dfj,{(p,(p^) (p (xi) ...(p {Xn) (pH^n+l) ■■■4'Hxn+m) &W iS [(p, (j)''] (1.35) 

What about writing this in terms of the fields (p' and (p' dagger? There are two ways of doing this. The first that comes 
to mind is to plug in the transformations of the fields, this gives, 

G{Xi,...,Xn,Xn+l,-;Xn+m) = ^ j D^{(p',(p'^) {cxp i {a (Xn+l) + ■■■ + a {Xn+m) ~ a (Xi) - ... - a (Xn))) 

(1.36) 

X (p' (xi) ...(p' [xn) (p'"^ {x„+i) ...(p'^ ixn+m) exp iS [(/>', (p>'^\ X exp - i / d^a; a ix) dy,f . 



However, as the fields are just dummy variables which get integrated over I should be able to write the generating 
functional, hy its definition , as, 

Z = j {<P', cp'^) exp iS [<p', <p'^] (1.37) 
and hence the Green's function, by its definition, as, 

G {Xi,.. ., Xn, Xn+l,. .., Xn+m) = ^ j D^{(p',(p'^) (p' (xi) ...(p' (Xn) (p'Hxn+l) ...(p'Hxn+m) CXp iS [(/)', (p'^] 

(1.38) 

The last two expressions for the Green's functions must be consistent. The last equation does not depend on the 
function a, 

= (1.39) 



Sa (y) 

so one condition for their consistency (there are undoubtedly others) is that G{xi, ...,Xn,Xn+i, ...,Xn+m) has no 
dependence on a. The previous statement is the key to understanding the origin of Ward identities and hence this 
course. Generally, in doing a transformation on some fields in a path integral one will introduce some parameters of 
the transformation into the expression, but given that the Green's function contains no such parameters by definition, it 
obviously can have no dependence on any such parameters, hence one can always apply the above condition. 

So what is the functional derivative in this particular instance? It can be useful when doing this to consider again 
the spatial coordinate as being just a label for the "true" variables that we are trying to differentiate with respect 
to. We are differentiating G with respect to the function a at a particular point in space yl To be idiot proof this 
means we get a delta function every time we differentiate an e*"'-^-' and differentiating e"'^"* ^ <^i^)9''ji^{x) gives 
^''i/i {x)\^^y e~^^ ''"^ a{x)d''j^{x) i^ggg jyvo types of term in the Green's function are the only ones that are acted on 
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by the differentiation as nothing else has any a' s in it. We are essentially differentiating the following, 

(exp i (a (x„+i) + ... + « {xn+m) - a (xi) - ... - a (a;„)) x exp - i / x a (x) {x)) 

= 7^ (sxP ^ (" i^n+i) + ... + a (a;„+™) - a (xi) - ... - a (a;„) - / d^x a (x) Q^j^ (x))) 

/ \ (1-40) 

= z (x„+i - t/) + ... + ^ (x„+„ - t/) - ^ (xi - y) - ... -S{xn-y)- d^j^ i.x)\x=y) 

X (exp i {a (x„+i) + ... + a {xn+m) - a (xi) - ... - a (x„)) x exp - i Jd^x a (x) (x)) 

= ^ J Df_,{4)',(f>'^) (exp j (q:(x„+i) + ... + a (x„+™) - a (xi) - ... - a (x„))) 

X z (^(5 (x„+i - y) + ... + (5 (x„+™ - y) - S {xi - y) - ... - S {xn - y) - (a;)!^^^) 

X (j)' (xi) ...0' (x„) cj)'^ (x„+i) ...(/-'■f (x„+„) exp iS [(j)', (j)''^] x exp - i / d^x a (x) (x) . 

(1.41) 

Note we haven't specified an a, this result is true for all a's! Let's see what happens with a (x) = 0, which clearly 
corresponds to not transforming the fields at all: 



5G 



Sa{y) 



= 

a=0 

= :! / 'A'^) i (a^n+i -y) + ... + 5 (xn+m - y) - S {xi - y) ~ ... 

-S (x„ - y) - a% J x (xi) ...<^' (x„) <^'t ...^'t exp iS [<j>', <j>'^ . 

(1.42) 

(s (x„+i - y) + ... - 5 (x„ - y) - d'f^y^j^ (y)j (xi) ...0' (x„) (x„+i) (x„+„) 



0)=0 (1.43) 



This last equation is a Ward identity. The terms 6{x — y) are known as contact terms they do not form an important 
part of our discussion and will play no further role. Consider the case where we have no fields before the exponent in 
the path integral i.e. m = n = in which case, 

(0 la^'i^l 0) = 0. (1.44) 

Considering this at the classical level one gets d'^j^ = i.e. one recovers a conserved current - Noether's theorem 
from the Ward identity. 

At this point a quick summary of the process might be useful. We took a simple scalar 0* theory with fields 
and (f)2 and rewrote it as a complex scalar field theory with fields <f) and (j^ this simphfied the form of the Lagrangian 
considerably. We then asked what happens to the Lagrangian when we make U {!) transformations of the fields. The 
result was that all the terms in the generating functional were invariant except for the kinetic part of the action which 
acquired a current term — J d^x a (x) Under the transformations the Green's functions were also not invariant 

and became messy under the gauge transformation a. However the definition of the Green's function contains no such 
terms and has the same form whatever the fields are that we are integrating over, the fields are "dummy" variables, 
hence for consistency we demand that the Green's function has no dependence on a when rewritten in terms of the 
transformed fields. Applying this condition to the expression for the Green's function in terms of the transformed fields 
yielded the Ward identity, which we shall see is all important in field theory. 

Classically the symmetry (Noether) currents were exactly conserved and the derivation of this depended on the 
equations of motion being satisfied. In the path integral formalism we integrate over all possible configurations not just 
the classical trajectory but we can expect an analogous result to that of the classical case, i.e. we expect the expectation 
value of the synnmetry current to be conserved as we know that the classical equations of motion hold as quantum 
averages. 

We will now attempt to generalize the formalism we just used in deriving our first Ward identity. First we define 
a set of fields (j)^ [x) 2 = 1, n. Next we define a local, infinitesimal unitary transformation Ay (x) = e {x)^ tf^ 
belonging to some group in which the action is synmietric, where are the group generators and e (x)^ are the 
associated infinitesimal parameters of the transformation. Finally we consider what correlation to study, we will be 
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general and say that we are looking at (O ({0*}) | O) where F is some function of the fields Under the 

transformation we have, 



(j)" (x) = (jf {x) + d<p' (x) 
representing a local change of variables in the path integral. It is evident that we have 



,iS[4>] 



(1.45) 



(1.46) 



since the result cannot depend on the name of the integration variable. 

We now make the all important assumption that the Jacobian to to change from D (0") to D ((/)*) is one. This 
assumption is non trivial since we are deaUng with compUcated functional integrals the meaning of which should 
really be carefully examined. It is often possible to prove that the Jacobian of the path integral measure is one and 
when possible we shall attempt to show it. In some cases the Jacobian is not one. When the Jacobian is not equal 
to one we generate anomalies which we will talk about in later lectures. In other words anomalies arise when a 
theory looks as if it possesses a certain symmetry due to the invariance of the action but it is not actually invariant on 
account of the non-invariance of the path integral measure. Until further notice we can forget about anomahes and take 
D (0") = (0') (...at least until lecture 4 - scale invariance). Continuing from the last equation we then have that, 



J D \f {{r}) e^^M - F {{^^}) e'Sl. 



= 



The variation is due to the functions e (x)"^ so we can rewrite the variation as. 



/dot/ 



d^.T 



Se {xY 



-e{xr =0 



note we integrate over all space-time to "add up" the variation everywhere and get the total. 



jD{f)Jd''xe{xflF{{^^}) ^^^^ 
Jd^xe{xfjD{c^^)\F{{ct>'}) 



SSI.P] , SF{{4^})' 



(1.47) 



(1.48) 



(1.49) 



As each of the functions e {x) are arbitrary then we require that. 



/ 



Dif)iFm)^ + ^^-ii^e^^l^ 



Se (x) 

5S [(p] 



Se {xy 



OVA 



Se (x) 



Se (x) 



(1.50) 



(1.51) 



For actions which are globally invariant under some symmetry the variation in the action produced by making the 
symmetry transformation local is the divergence of the symmetry current, for unitary transformations this is d^j^ oc 

Now we shall consider a "higher" synnmetry, an SO (3) global synmietry to illustrate more features of the Green's 
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functions. Consider the following triplet of fields and their transformation by R€ SO (3): 




i?(/>. (1.52) 



From this representation we can easily construct a globally SO (3) invariant Lagrangian: 

T T T 

Kinetic Mass Potential 

term term 



(1.53) 



In our path integrals the measure of the integration is clearly invariant in exactly the same way as it was for the fields 
(/>! and (j)2 in the example at the end of lecture 1. 

jDct = n./_rd'Ai(^)/_rd<^2(a=)/+^d<^3W 

= n. /.r #i (^) /-r #2 (^) /.r #3 (^) >< oet (i?) 

= 

Note that if this was a local transformation we still have no problem as we would get the same product of ones 
(Det {R {x) ))'s over all space-time points x. However if we were dealing with a local O (3) transformation instead then 
Det (i?) could be —1 or +1 and we end up with an ill-defined product of +l's and — I's over all space-time points, this 
means that the generating functional and Green's functions would be badly defined in terms of the transformed fields. 
It is perhaps worth noting that we stick to special (Det = +1) groups in our course for this reason. 

It is a general result that only combinations of the fields for the Green's functions arguments which have a compo- 
nent that is a singlet representation of the transformation group are non- vanishing! For instance consider the following 
vector of Green's functions, 

G{xi) = ^jD^c^,ixi)cxpiS\£ 
^ G{xi) = |/£l£i?H0^(a;i)expiS'[£] . 

As the action and integration measure are SO (3) invariant then the following integrations/Green's functions should all 
be equal, we can show this if need be with a simple change of variables (no Physics just Maths), 

G {xi) = ^Jd^ ^1 (^i) cwiS[^=^J D£ <A'2 {xi) sxpiS[^=^J D£ exp iS . (1.56) 

This gives, 

G{xi) = {Ru + Ri2 + Ria) ^ -^i (^i) e'^P {Ru + R12 + R13) G{xi). (1.57) 

As Rij G SO (3) is arbitrary we have that Ru + R12 + Ris is not generally equal to one thus for consistency in the 
definition of the Green's functions it must be the case that G (xi) = 0. This result can also be obtained by imposing the 
(usual) constraint that the Green's function (by definition) cannot depend on any parameter(s) of the transformation. 



R 



^ sinpsimp —coscpsinip —cosip 

cos(f)sin6 — cos6cosij)sin(j) cos6cos(f)costp + sinOsimp —cosOsimp \ (1.58) 
\ cosOcoscp + cosipsin0sin(/) —coscpcosipsinO + cosOsincp sinOsinip 
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5G {xr) _ 6G {xi) _ 5G jxi) 

se ~ Sip ~ 5(j) 



= 0. 



(1.59) 



Note SO (3) transformations only have three parameters ( = number of generators) not as many as the number of 
elements in R so they can only give at most three conditions by this method. 

The property of non-singlet combinations of fields having vanishing Green's functions was also seen in the C/ (1) 
case just considered where we saw that only strings of fields (before the exponent in the path integral) containing an 
equal number of ^'s and ^^'s were non-zero. So in the case of SO (3) this means we are interested in things like, 



When dealing with Green's functions in this way it is important to realize that just because the string of fields isn't 
invariant under the symmetry transformation does not rule out the fact that it may have a singlet component (i.e. it 
does not mean it is equal to zero) it could have a singlet component inside it somewhere. One must project out the 
irreducible representations of the symmetry group from the string of fields preceding the exponent in the path integral 
and then look for the singlet component. 



G(xi,X2) = J D4> (f^ {xi) 4> (x2) exp iS \4i\ 

= f D(j)' (j)''^ {xi) R^R<t>' {X2) exp iS \ 
= jD£^{xi) .£ {X2) exp iS . 




(1.60) 
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1.3 Gauge Symmetries. 

In this lecture we generalize our knowledge of local U {1) symmetries to more complicated symmetry groups. Consider 
the following general unitary symmetry transformation , 

(t>a ix) ^ ix) = Uabh {X) (1.61) 

UlUkj^S^j (1.62) 

where a repeated roman index implies a sum over it. We shall restrict ourselves to special unitary transformation i.e. 
Det {U) = +1, this will avoid the nasty problem of getting ill-defined products of — I's and +l's in transforming the 
integration measure. 

The kinetic term in the action for our system of fields transforms as follows, 

If the transformation is a global one then the U^^ and Uac terms can be brought in front of the derivatives to give She- 
lf we take the potential to be once again a function of (pl^(f)a as before then both it and the mass term m?(l)\(j)a wiU be 
invariant whether the transformation is global or not: 

<^I<^a = <l^'Mai.^)Uac{x)<l>', 

= <l>'b^StcCb', (1.64) 

Next we wiU study the effect of making this general (special) unitary synometry a local symmetry. To begin with we 
define the symmetry transformation in more detail. Writing a general transformation in terms of its parameters ua {x) 
and its generators tA we have, 

U (x) =expiJ2aA{x)t^. (1.65) 

A 

The generators are defined by the following relations, 

jABC gj-g numbers known as the structure constants of the group. Now we rewrite the kinetic part of our action in 

terms of the new variables cj)'^ = U~f^ {x) 0^. For brevity I will drop the indices on the fields and transformation matrix 
so from now on (t>\U~^ {x) is (j)^U~^ {x), Uab (x) 4>b is U {x) (j) etc... 

5fe [C/t<A',<^'tc7] = /d^a;<-a^(<^'tc7)a,(C/t0') 

= / d^x 7?^'' ((9^<i6't) U + cl>'^d^U) [Wd^cj)' + (9,[/t) cj)') 

= /d^xry^-'a^^'ta^../)' (1.67) 

+ / d^x 77^"^ { (a^(/.'t) {ud,u^) (/.' + (/.'t t/t5,0'} 

+ /d^xr?'^''</.'t(5^C/)(a,C/t)(/.' 

Now we need to have a think about terms Uke 9^f/ {x). We omit the Y^a^ repeated indices to represent sums. 

d^U {x) = Q^exp iuA {x) 

= d^{l + iaA{x)t^-l{aA{x)t''){aB{x)t'')+...) (1.68) 
= i{d^aA{x))t^ + 0{a'') 
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So for infinitesimal transformations we can safely drop terms of order ua {x)^ and above. Back to the action, the 
second term, 

J d^x T]^'' { {d^(t>'^) (Ud^U^) ct>' + 4>'^ [d^U) [7ta,0'} (1.69) 
neglecting terms of order a a [xf and above becomes, 

-i J d^x <^ {(a^</.'t) {d^aA {x))tW - <^'^ {d^aA {x))t^ {d^(f>')} . (1.70) 
If we exploit the symmetry of rj^i, we get, 

- J d^'x <^ {d^aA {x)) i{{d,cl>'^) t^c/,' - ct>'h^ id,cl>')} . (1.71) 

We now define our current, 

= i t^</.' - {d,^') , (1.72) 

which makes the middle term in the action equal to, 

- j d^x {d''aA{x))j^ = - j A^'xd'' {aA{x)j^)+ j d"" x aA (x) d'^j^ . (1.73) 

On the right hand side of equation 1 .73 is a 4-divergence which we can express as a surface integral. If all the fields ^' 
and (p'^ inside vanish at infinity (the surface would be integrating over) then so does and the surface integral is 
zero and, 

- J d^x {d'-aA {x}) 3^ = j d^x aA {x) d'' (1.74) 
The kinetic part of the action is then, 

= Sk [</-', + /d^x aA (x) d"j^ + /d^x n^-^'^ (d^U) {d,W) 0'. 

As we are neglecting terms of order aA (x)^ and above the last term in 1.75 is effectively zero since its first non-zero 
contribution is of order aA (x)"^- This gives us that the entire action is invariant up to a current term which arises from 
the variance of the kinetic term. 

= Sk[cl^',4>'^]+Jd°xaA{x)d-j^ 

From this one obtains the Ward identities in the usual way by functionally differentiating with respect to the param- 
eters of the transformation. Recall that we must have some singlet representation of the fields in the Green's function 
for it to be non-vanishing. This requires that we have an equal number of ^ and (j)^ fields (this is essentially stating 
global charge conservation). 



G{xi,..,Xn,Xn+l,--,Xn+m) = ^ J {(f), (f)^) (p^ {xi) (j) {Xi) ...(j)^ {Xn) <j> {Xn) exp iS 

= ^JD„ (f/^) (/)'t e-i«^(-"+i)*-*e*"^("i)*^</.' (xi) ... 

X 0't {Xn+m) e-'"-^(^"+'")*''e'"'^(^'*)*''?i' {Xn) CXp iS [(p, ^t] 

X exp / d'^x aA (x) d^j^ 



(1.77) 



As usual, by definition of the Green's function we have also, 

G (Xi, .., Xn, Xn+l, --, Xn+m) = ^ [ Dfj, {(j), (j)^) (p^ (x„+i) (j) (xi) ...(/)^ {Xn+m) 4> (Xn) exp iS [(p, (p^ (1.78) 
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which has no dependence on any of the (a;)'s. Functional differentiation with respect to the parameters a. a {x) and 
setting (a;) = V A gives the famiUar looking form for the Ward identities: 

(O I ((5 {xn+i -y)t^ + ... + S {xn+m - v) - S (x, - y) - ... - 5 (x„ ~ y) - d'(y^3^ (y)) 

X 4>^ {Xn+l)4'{xi) ...(j)^ {Xn+m)(t>{Xn)\^) = 

(1.79) 

An obvious question to ask now is what do these currents couple to? To answer this we invoke the gauge principle 
and introduce new gauge fields to make the Lagrangian completely invariant. Given that our fields (j) transform as 

(j) ^ (f)' = U {x) <f) and 0^ <f)'^ = <f)^ {x) this means, 



d^cf>' = {d,u (x)) <j> + u (x) {d,<j>) 
a,<^'t = {^^(t>^)W{x) + c|>^^,u^x)) 



(1.80) 



which we can rewrite as, 



(1.81) 



d,^' = u{x){uHx){d,u{x))d^ + d,(i)) 
a,(^'t = {{d^(i>^)+(j,^{d^W{x))u{x))W{x) 

Sandwiching these two terms together gives a kinetic term which isn't gauge invariant. 

a,</,'ta-^' = ((a^^t) + 0t ^Q^u^ (a-)) u {x)) [u^ {x) [d^'u [x)) <j) + d"^) (i.si) 

To make this gauge invariant we need to introduce a gauge field which transforms as A'^ = U (x) ^^^t/^ (x) — 
{dfj^U {x)) W (x). This gauge field is added to the normal derivative to form the covariant derivative -D^u = dfj^ + A^ 
so-called because it transforms like the fields do: 

D'^<f>'{x) = d^<f>'ix)+A'^ct>'ix) 

= {d,,U{x))cbix) + Uix)d^<f>ix) + Uix)A^<l>ix)-{d^U{x))<l>{x) 
= Uix)id^ + A^)c^{x) 
= U{x)D^(j){x) 

Note that in future we will not write down the explicit space-time dependence of U,U ~ U (x) unless otherwise stated. 
If we now replace the derivative in the kinetic part of our action by the covariant derivative we will find that we 
have something which is now gauge invariant i.e. replace 

v'"'d^(l}^d^<l}^V^''{D^<f>)^ D^<P. (1.84) 

It is worth making a few points about the field A^ that we have introduced. Firstly A^ is actually an anti-Hermitian 
matrix. The anti-Hermitian property is easy to prove, imagine = 0, this is called a pure gauge configuration as 
gauge transformed versions of the field consist only of the gauge transformation's inhomogeneous term, then 

a; = UA^W-{d^U)W 
=> A'^ = U{d^W)-d^{uW) (1.85) 

= u{d^W)-d^{i) 
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Also for global transformations {d^U — 0) is a matrix living in the adjoint representation of the transformation 
group i.e. in this case we have A'^ = UA^W. In hght of this fact it is possible to write the matrix in terms of 
the generators of the group (x) = A/^^a (x) t^- Note that in the case of the Abelian gauge field the covariant 
derivative appears to contradict the definition above, in massless QED we have £ = ■tpi /D-ip = {i /)+ /A) tp 
i.e. Dn = d/j, — iAfj, which implies that the gauge field transforms as — > UA^W — i {dfj,U) W, this is just 
the usual — > + (x) for U = e"'"^^^. This pecuUarity is an annoying reUc in the general Uterature and 
arises because as the Lagrangian there the kinetic term has an i attached to it but the gauge field does not. With this 
definition the field in QED is Hermitian, this is easy to see as one just has an extra factor of i floating around in 
the above calculation where we proved it was anti-Hermitian. Henceforth we use the ip [i /)+ /A) ijj form of the QED 
Lagrangian. A better definition of the QED Lagrangian would be ipi ( /)+ /A) ip in which case A^ does transform 
as A'lj^ = UA^W — {d^U) W and we have a better analogy with the Lagrangians and gauge transformations of 
non-Abehan gauge theories. 

A^ is sometimes referred to as a connection which means it does for gauge theory what the Christoffel symbols, 
r^^, do for General Relativity. InGeneralRelativity the covariant derivative of a vector is D,yV^ = d^V^+T'^^V^ 
and in our gauge theory we have Z),y(/)" = 9,^(6" This is perhaps interesting from the point of view that General 
Relativity is a gauge theory locally invariant under transformations in the Poincare group. It is also worth noting that 
doesn't transform as a tensor except for global Poincare transformations (it has an inhomogeneous term in its 
transformation law), 

,^ _ dx^dx^dx^ d'^x" dx'" 

~ fe" dx'^ dx't" dx'^dx't' fe" ■ ^ ' ^ 

For global Poincare symmetry we have ^73- is the Lorentz transformation matrix plus some matrix of translations i.e. 
dx'i^dx'-y ~ ^ /3) 7 and hence r^_^ would transform as a tensor, the inhomogeneous term is identically zero. Now 
if r^;^ transforms as a tensor by the equivalence principal it is possible to transform to a frame where it would be zero 
at a given point and hence it would be zero everywhere i.e. T^^ wouldn't exist - in special relativity F^^ is effectively 
transforming as a tensor, it doesn't exist. This is analogous to our gauge field A^ it also has an inhomogeneous term in 
its transformation law — {d^JJ) . 

The only thing left to do now is turn our (x) field into something physical. To do that we have to give it a kinetic 
term in the action. The kinetic term should be gauge invariant. The following relation is essential to constructing these 
terms. 

= {d^ + A'^){UD,cb) 

= {d^ + UA^W -id^U)W){UD,(b) (1.87) 
= (d^U) D^cj) + Ud^D^cP + UA^D,<t) - {d^U) D,<j) 
= UD^D^cj) 

=^[D'^D',](j>' = U[D^,D,]cj> (1.88) 
We denote F^^ = [D^, D^] = dvA^ — d^A^ + [A^, A^]. From the above we see that, 

F^.<P^F'^^<P=UF^^cP = UF^^U^U<p (1.89) 

from which it is clear that F^i, transforms as, 

F^.^F'^, = UF^^Ul (1.90) 

To make the A^ field physical we then take its gauge invariant, Lorentz invariant kinetic term to be 

Tr [F^.F"-] ^ Tr = Tr [UF^^WUF^^'^W] 

= Tr [UFi,^F^"'U''] . 
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It is a result that the trace of a product of matrices is independent of the ordering of the product hence, 

= TrlF^^F^^-^] 

This means Tr [Fni^F^^"] is indeed gauge invariant, 

Tr [F^^i^'A"^] = Tr[UF^.F^-^W] ^^^^^ 



Tr [F^.F^"' 



Chapter 2 

Conformal Symmetry. 



2.1 The Conformal Group. 

This lecture is about space-time symmetries. We start off with the usual scalar field action, 

S [(!>] = j dPx ^d'^^d^^ - m'^cf' - y (0) . (2.1) 

In natural units {h = 1) we have that S [</>] is dimensionless. For the dimensions of the terms in the action to be 
consistent we then require that m. have dimensions of length^^ as 9^ has dimensions of lengthr^. This action is 
invariant for transformations of the Poincare group i.e. it is invariant under spatial translations, rotations and Lorentz 
boosts (we are integrating over all of space-time). As a result of these symmetries we already have some important 
consequences. Consider the effect of translational invariance on the Green's functions for instance: 

G{xi,...,x„) = ^ J D(f) (f){xi) ...(j){xn) expiS[(j)] 

(j) (Xf,) = 4>' {x^ + a^) (2.2) 
S[<t>\ = SW] 

The measure is also invariant, consider again the spatial coordinate x as the index labeling the true variable in the path 
integral i.e. field. The translation amounts to re-labeling this index (by a constant shift), the path integral measure is 
the product of the integrals over all possible values of the field at each (and every) space-time point so we do not expect 
the measure to change, roughly speaking, 

nxd</'(^) = n.-a#(^-a) 

= n.-a#'(a:) (2.3) 

= U^dcP'ix) 

=^ G{xi,...,Xn) = ^ J D(j)' (j)' (xi+a) ...(l){xn + a) expiS[^'] ^^^^ 
= G{xi+a,...,Xn+a) 

So for an action and an integration measure which are invariant under spatial translations one has that the Green's func- 
tions are also translation invariant and hence they can only depend on where their arguments are relative to each other 
i.e. Differences in the positions of their external points. Similarly one can deduce from Lorentz/rotational invariance 
that the Green's functions can only be functions of things like (a;i — X2)^ {xi — 2:2)^, {xi — x^Y {^i ~ ^^c... 

In actions with m = we have an additional "special" symmetry to consider for particular choices of V {(p). We 
start this study by considering the dimensions of the various pieces of the action. Let the fields (j) have dimensions 
of length^ this means that for the kinetic part of the action we have overall dimensions length'^^~'^^^ , we have a 
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product of two fields giving length?^, two partial derivatives, length~'^, and an J d^x giving length^. At the start 
we said the action was dimensionless so this means A = The potential term should also be dimensionless. Say 
V {(f) ~ g4>^ with g a dimensionless coupling constant, this has dimensions length^^"^^^^ but when we add in the 
effect of the integral / dPx the potential part of the action actually has dimensions length^ which of course 

must actually be dimensionless i.e. P = 2D/ (£) — 2). So in £) = 4 space-time dimensions we have V (</>) = gcp"^, in 
D = 3 space-time dimensions we have V {(j)) = gcj)^. 

The "special symmetry" which we alluded to earlier is a symmetry of the action under scale transformations 
x'^ = Aa;'*. Under such transformations the dimensionless action that we deduced above transforms as follows. 



S[<j)] = Jd^x ^d^'(l){x)d^(l){x) - g(j){x)-^ 

= A-^ / d^x' A^is'^c/, (^) a;^ (5^) - g^ ■ 



(2.5) 



where &, 



/ _ d 



jii Xdxf • 



2D 
D-2 



If we now define that the fields transform under scale transformations as (j)' {x') = A~ * 2 ' ^ t^eu have the 
action is invariant under scale transformations: 



S[(^] = Jd''x^d^'c^{x)d^<P{x)-gct>{x)^ 



^3^'c^{x)d^(P{x)-g<t>{x)^ 
= Sm (2.7) 
= / d^x (x) 5^<^' {x) - gct>' (x)^ 

For these scale transformations the measure is not invariant! 

= n..A^/d<^'K) (2.8) 

= (n.,A^)/w. 

However at the front of the definition of the Green's functions we have a factor ^. The generating functional Z acquires 
exactly the same factor when the measure is transformed [Yix ) in exactly the same way and so the two factors 



in the numerator and denominator cancel. As a result of this scale invariance the form of the 2-point function is fixed! 



G{xi-X2) = i/D(/.A^(/)'(A,xi)A^0'(Ax2)e^^K(^')] 

= A^-2i/i^,/),/)'(Axi)<^'(Ax2)e'^K(^')] (2.9) 
= A^-2G(A(xi-X2)) 



This last equation in fact means that the 2-point function must be of the form, 

G{x,-X2)^- ^—jj^. (2.10) 

\Xi - X2\ 

This result is true in Minkowski and Euclidean space. 

So far we have considered scale transformations and Poincare transformations separately. The combination of the 
Poincare transformations and scale transformations are known as the conformal transformations. Note that Poincare 

and scale transformations preserve the angles between things. Here when we talk of angles we are talking about angles 
in whatever D dimensional space we are in i.e. if we are in Minkowski space we mean the angle between things in 
3 + 1 space-time dimensions not 3 spatial dimensions. Perhaps a better, more formal definition of a conformal (angle 
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preserving) transformation is that it is a transformation where the ratio, 

— ^ (2 11) 

\dx\ \dy\ 

is unchanged. Clearly rotations, translations, Lorentz boosts and scalings all satisfy this requirement. 

Under the general (infinitesimal) coordinate transformation x ^ x' = x + e{x) the metric tensor undergoes the 
following transformation to first order in e {x), 

— yap dxt^ (Jk" 

= 9o.p{5'^ + d^e"{x)){5P + d,e0{x)) (2.12) 
= 9u.v + gixffdi^e^ (x) + Qai^d^e" {x) 

If it is the case that the transformation in the line above is of the form g^j^^, g'^^ = g^^ x (1 + A (x)) then the metric 
is just being multiplied by some number (at each space-time point x) and the transformation is a local conformal 
transformation i.e. it is (locally) angle preserving. We would Uke to find functions e {x) that have this property. Let 
us take the case that our initial metric is flat g^j^v = r]^v in which case we want di^e^ {x) + d^ev {x) = A {x) rji^v 
Contracting the last equation gives 2d^e^ {x) = DA (x) where D is the number of space-time dimensions we are 
considering. Substituting in for A (x) we have, 

2 

^ d^Cf, (x) + d^e^ (x) = -5'=e^ {x) r]^^. (2.13) 



differentiate this with respect to 9'* to get. 



d^d^e, (x) = ( I - 1 ) 5, {d'^e^ {x)) , (2.14) 



and differentiate again with respect to d", 

2 



^ 2j d'^d^ {d^'e, {x)) = 0. (2.15) 

This equation has two solutions. Either D = 1 or d^d^ {d^Cv (x)) = 0. To solve for (x) consider 2.13 and 2.14. 
Acting on 2.13 with d^d\ we get, 

d^dxd^e^ (x) + d^^dxd^e, (x) = ^d^d^d'^e^ (x) r,^^ (2.16) 

but the equation which we are trying to solve tells us that d^dxd'^e^ (x) = so 



d^dxd,e^{x) + d^dxd^e^ix) = 

d^d^dxe^ix) = -d^d^dxe,{x) 



(2.17) 



If we now act on 2.14 with 9^ we find. 



d^d^dxe, {x) =(^^-l^ d^d, {d^ex [x)) . (2.18) 

Looking at the above equation we see that the right hand side is symmetric in /i and v so the left hand side must be 
also however we have just shown that d,yd^dxe^ (x) = —df^d'^dxf-v ix) so the only way that we can have these two 
equations consistent is if = 2 or d^d^dx^v [x) = and hence {d^ex (x)) = 0. 

We shaU not discuss the case where D = 2 but try and solve d^di, (d^ex {x)) = 0. With a littie thought one sees 
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that the general solution to this equation is d^e\ {x) = A + B^x\ which means {x) must have the general solution, 

{x) = + hni^x" + Cn^xx^x^. (2.19) 

For a non vanishing x"x^ term we must have c^vx symmetric in u and A. This is trivial to show, first use the fact that 
V and A are dunmiy variables, summed over, so we can call them whatever we Uke, we can call u A and call A u 

=^ CjxvXX X^ = CfiXi/X'^X 

Now we have a form for the e^^ (-7^) 8 where the e^^ (x)'s are performing a rescahng of the metric. If we plug our 
general solution for {x) back into our original equation for {x) i.e. 

2 

d^e^ {x) + d^e^ {x) = -a'^e^ (a;) 7?^^ (2.21) 
we can get constraints on the coefficients and c^vx- First consider 9i,e^ {x), 

dv^n [x) = (o^ + b^pxP + c^pffxPxi^) 

= {hp^pXp + Cp^piixPx^) 

= b^p5P + c^pii5PxP + c^p0xP5i (2.22) 

— b^ij {^fii^p xP 

— bp,!/ ~t~ '^Cpi/pXP . 



Substituting this into 2.18 gives. 



2 a ** a 

^/ii^ + "^CpupX + hup + 2c,jppX = -^bxTjpv + — c xpX Vni^- (2.23) 



Comparing the coefficients of xP on either side of the equation one gets the following constraints on the forms of bpv 
and Cp^x, 



bpu ~\~ bvp — j^TJ^iv^x 

2 A 
C-nvp + Cvup = 'DVpf'^ Xp- 

Let us write bp^ in the form, 

= Mp^ + ^Vp-i'- 
Substituting this into bpi, + b^p = -j^Vtiv^x ^® have, 

Mp^ + %r]p^ + M^p + %r]^p = %Vi.^^b^ 
Mp^ = -M^p 

i.e. we can write bpi, as the sum of an antisymmetric matrix Mpi, and a symmetric matrix rjpvA, 

bpv + bt,p = Mp^ + r]pi,A + M^p + r]„pA 

= iM + M^)^^ + 2r^p^A 

= 2r]p^A 

= '^VpiyT)bx- 



(2.24) 



(2.25) 



(2.26) 



(2.27) 
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T T 

Lorentz Scale (2.28) 

Transformations Transformations, 
[incl rotations) 

Turning to the constraint on c^i^a we have, 

c^vp + Cvnp — '^Vfj.i'^ Xp- (2.29) 

We now wish to constrain c^i^x with an analysis similar to b^,^. We can do this by remembering that a non- vanishing 
x^x"^ term in (x) requires that c^i^a is symmetric in i/ and A. 

2 A 

: , _ (2.30) 

= %Vfiv(^\p ~ T)'>li^pc\n + -hVpiiC^Xj, — C^„p 

=^ c^i/p — j-^Vp.v'^ Xp jj'Hvpc Xp + jjVpnC xu 
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2.2 Conformal Symmetry and the Dilatation Current. 



In this lecture we round up our study of scale transformations. From last lecture we had that for have a general 
conformal transformation, 

T T T T 

Translations Scaling Rotations/ Special Conformal 

Boosts Transformations 

(2.31) 

If our theory is symmetric under such transformations then associated with the symmetry should be Ward identities 
and conserved currents as in the previously studied theories. We expect to find Ward identities for both rigid rescaling 
A ^ as well as the special conformal transformations. 

In the last lecture we saw that the following action was invariant under rescaling a; — > a;' = Aa; (A is a constant), 



(2.32) 



with the fields transforming as (p' {x') = A 0-2 0(a;). We shall now make the following infinitesimal change in 
coordinates (we neglect terms less than of order e (x)^), 

x'^ = {l + e{x))x^ (2.33) 

(j)'{x') = {l + e{x)y^ (f>{x). (2.34) 
The Jacobian of this transformation is needed for the integral measure d^x in the action, it is given by. 



J (a;, x') 



dx'" 



S^^{l + e{x))+x^'^^e{x)\ 
1 + e (a;) + x^die (x) x^d2e {x) x^d^e {x) 

x'^die (x) 1 + e (x) + x^d2e {x) 

x^dic (x) 



l + e{x)+ x'^doe {x) 



(2.35) 



With a little thought one can see that all of the off diagonal terms in the determinant above will only contribute terms 
of order e (x)^ and above so we can forget about them at all points in evaluating the determinant making it the product 
of the diagonal terms. To first order in e (a;) the determinant is. 



J {x, a;') = 1 + De {x) + x^d^e (x) +o(e {xf^ . 
Jd''x=J d^x' jj^^ ^ J d'^x' (1 - De {x) - x^d^e (x)) 



(2.36) 
(2.37) 
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We attempt now to evaluate d^4> in the new variables so as to ultimately work out the transformed action. 

« (1 + e (x)) + xP'-0) ^ ((1 + {x)) ^' {x')) 

= {5^^ (1 + e (X)) + + o^, (,) £gLl~ 

+ (1 + e {x)) + xpMeI) d^^> (^') M£) 
= ^^(l + e(.))(^ + ^e(.)£|M)+,P^M|) (2.38) 

~ dx'" ^ 2 ^y-^) ax'" ^ " V"^/ Sx'" ^ 9x'P ax" 

_ d^'{x') D^x , X 9'/''(^0 , p90V) M£) + ^D^rti' rr'l ^ + O (f^N 

~ ax'" ^ ^ 2 j t i^x; g^,^ -I- X -r 2 Y J ax'" ) 

= i^ + hix))%^+x^%^'^ + ^<p'ixr^ + 0{e^) 

For ease of notation we denote differentiation with respect to the transformed coordinates x' gpjr as d'^ i.e. 

d^4> {x) ={l + ^e {x)\ d'^4>' {x') + x^cl^' {x') d^e (x) + ^-^ct>' {x') d'^e (x) + O (e^) . (2.39) 



This means the kinetic term in the Lagrangian is, 

^d>^ct> (x) 5^</> (x) = \ {d'^^' ix') + f e (x) 5;^' {x') + x^d'^^' (x') d^e (x) + {x') d'^e (x)) 



X [d'f'cj)' {x') + f e (x) d'i'cl)' {x') + xPd'pCj)' {x') d^^e [x) + ^(f>' {x') d'l^e [x)) 



2' 

\^'^'(|)' {x') d'^cj)' (x') + \De {x) d'f'cj)' (x') d'^(f>' (x') + xPO'^c/)' (x') d'^(l>' (x') Q^e (x) 
+^0' (x') (x) d'^(f>' (x') + O 



= i (x') (x')) (l + De (x)) + (S;^' (x')) (x^a;^' (x') d>^e (x)) 

+ {d',^'{x')){^cl,'{x')d"^e{x)). 

(2.40) 

In the action all terms in the Lagrangian will be multiplied by the Jacobian j^J"^^,^ = 1 — De (x) — x^'d^e (x) + O (e^) . 
So we have even more terms of order and above to find and then drop from the kinetic part. 

/ A^xld^'cj) (x) d^(j> (x) = / d^x' (1 - De (x) - x^dxe (x)) {\ {d'^cj)' (x') dlcj)' (x')) (1 + De (x)) 

+ {d^ct)' (x')) {xPd'p(P' (x') ^^'e (x) + ^</.' (x') S'^'e (x))) 
= / d^x' (I (x') a;^' (.x')) (1 + De (x)) + (5;,/.' (x')) {xPd'p^' (x') (x) 

+ ^(f>' (x') a'^'e (x))) - i (x') (.x')) {De (x) + x^^ae (x)) 

= / d^x' i (x') (x')) + (x')) {xPd'^^' (x') a^e (x)) 

+ (x')) (^</'' (x') &^e (x)) - ix^^Ae (x) (S'^c/.' (x') (x')) 
= / d^x' i (a'"./)' (x') a;<^' (x')) (l - a;^aAe (x)) 

+ (dlcj)' (x')) (a;^^;./)' (x') ^^'e (x) + ^</.' (x') d't'e (x)) 



(2.41) 



For the potential part of the action we have the following (dropping terms O (e^) and above): 



2D 
D — 2 



r/d^X(/)(x)^ = -5/d^x' (l-r>e(x)-x^aAe(x)) (^(l + e(x))^?!)'(x')j 

= -5/d^x' (l-£>e(x)-x^5Ae(x))(l + e(x))^0'(x')^ . (2.42) 

= -gjd^x' {l-De{x)-x^dxe{x)){l + De{x))(j)'{x')'^ 

= -gj6^x' [1 - x^dxe (x)) (f>' {x')^ 
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So the effect of the transformation / change of variables on the action as a whole is, 



J d^a; C{ct>)= J d^x' C (<^') (l - x^S^e (x)) + (5;<^' (x')) (^x^d'^c/)' {x') S'^e (x) + ^^^<^' {x') 5"^e (x) 

(2 

^ 5 [<^] = 5 [</>'] + y" d^x' - L (</)') x^S^e (x) + (a;<^' (x')) (^s^a^^' (x') a^e (x) + ^^^4>' i^') 5"^e {x)^ . 



(2.43) 
(2.44) 



To first order in e (x) we have e{x) = e {x'}. This is shown by Taylor expanding e (a;'): 

e{x'i^) = e{x^' + e{x'')x'')=e{x^')+ ^^P- e {xi") x" + ... 



■: (xi^) + o(e {x^'f^ . 



yU.=Xl^ (-2.45) 



Also, 



dle{x') = §^dxe{x) 

= l^{l + e{x))-\'^)dxeix) 

{^{l-e{x))x'>;)dxe{x) . (2.46) 
d\e (x) 



dx'" dx'" 



If we now use these identities on the second term in the transformed version of the action, equation 2.44, 

j d^x' - C {<!>') x^d^e (x) + [d'^cj)' {x')) (^xP&pCj)' {x') d'^e (x) + ^^<f>' {x') d'^e (x)^ , (2.47) 

we get, 

fd^x' - L {(j)') x^'^'^e {x') + {d'^cj)' [x')) {xPd'pCj)' {x') a'/^e {x') + ^cj)' [x') fff^e {x')) 
= J d^x' - C {(j)') x^^'^'e {x') + {d'^(j>' {x')) {xPd'p(f>' {x') + (x')) ^'^'e (x') (2.48) 

= / d^x' {-£ {,4>') X, + {d'^4>' (x')) {xP&^cf' {x') + i^<i>' {x'))) d'^e (x') . 

If we use the product rule for differentiating / integration by parts we can write this as. 



= / d^x' d'P ( (-£ (0') X, + {d'^<t>' (x')) (x^a;^' (x') + l^cf' {x'))) e (x')) 
- /d^x' e (x') 9'^ (-£ {^')x^ + {&^4>' {x')) (x^a;</.' (x') + (x'))) • 



(2.49) 



As the first term is a D-divergence of some function of the field then if we make the usual assumption that all fields 
and their first derivatives etc vanish at infinity then the first term goes to zero (we write the volume integral of the 
D-divergence over all space-time as a surface integral, the surface being at infinity). So the transformed version of the 
action is, 

5 [0] = 5 [<!>'] - j d^x' e (x') d'P (^-L {ct>') x^ + (a;<^' (x')) (x%^' (x') + (x')) ) • (2.50) 

As x^ = (1 + e (x)) x^ and we have an e (x') at the beginning of the whole integrand, neglecting terms of order and 
above, we can just swap x^ and x^ for x^ and x'P; 

S[^] = S [<!>'] + J d^x' e (xO a'" (</)') x; - (a;^' (x')) (x'^d'^cf,' (x') + ^<^' (x'))) (2.51) 

^S[(t>]=S [</.'] + J d^x' e (xO a% (x) (2.52) 
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where we have defined the dilatation current, 




(2.53) 



The dilatation current is quite novel on account of it having x'^'s in it. The dilatation current is a common feature of 



field theories incorporating gravitational effects, it couples to a scalar field, the dilaton. 

The Ward identities are obtained in the exact same way as before, we consider a general Green's function, make a 
change of variables (do the transformation) and then differentiate with respect to the parameter of the transformation. 
The measure of the path integral is invariant for exactly the same reason that it was in the last lecture (we pick up 



Up to now our discussion of our scale invariant scalar field theory has been at a pseudo classical level. If we are to 
think about the theory as a quantum field theory we will have to think about UV divergences. Consider for instance the 
first order correction to the two-point function in theory, this is a loop on top of the free particle propagator. To make 
our theory meaningful we have to be able to deal with the infinities produced by integrating over loop momenta, we 
have to renormalize it, which means introducing a dimensionful parameter, a cut off, into the otherwise dimensionless 
theory. Introducing dimensionful parameters into our theory breaks scale invariance. An example of this can be seen 
by trying to introducing mass terms into the theory. Consider the action of lecture 4 with a;'^ = Aa;'' for constant A, 



factors of (l — -^^j^e (a;)) in the measure which cancel with those picked up in the measure of ;^). 



S[(j)] = J d^x (x) df,(f> (x) - m'^cj)'^ (x) 




(2.54) 




(2.55) 



= Jd^x' ^d'''(b'{x')d'^(l)' {x')- - 



(2.56) 



Chapter 3 

Principles of Gauge Field Theory 
Quantization. 



3.1 Faddeev-Popov Gauge Fixing and Ghosts. 

This lecture deals with how to quantize gauge theories, specifically we describe how to get around the problem of 
integrating over too many/artificial (gauge) degrees of freedom in the path integral by the Faddeev Popov gauge fixing 
procedure. 

We will illustrate Faddeev Popov gauge fixing on QED though the method generalizes to the other gauge theories. 
The problem is that in the path integral formalism transition probabilities are obtained by integrating over all possible 
physical paths, but if we naively write the generating functional of QED as an integral over the gauge field A^^ we are 
integrating over more than just the number of physical paths because all fields related by a gauge transformation 
represent the same physical configuration (they lie on the same gauge orbit i.e. they are related by a gauge transforma- 
tion). They give rise to the same E and B etc fields that we observe. There are actually an infinite number of gauges to 
choose from which sure enough means that we over count the number of paths by a factor of infinity, ideally we want 
to isolate and divide out this "infinite group volume factor". This means we must somehow constrain our integral over 
the gauge field such that we are only integrating over all possible physically inequivalent paths. 

To solve this problem it is tempting to insert a condition into the path integral like say S (Lorentz gauge 

condition) or 6 (v.Aj (Coulomb gauge condition) which would set the contributions from all other gauge transformed 
versions of to zero i.e. we want something to fix the gauge (we only want to cut across a gauge orbit once when 
integrating through the function space of A^). Ideally we should do general derivations and so we should consider 
inserting a general gauge condition F {A^) = with F an arbitrary function of the fields (and their derivatives etc). 
However we cannot trivially insert such delta functions as this changes the measure of the integration DA^j^. Say we 
had an ordinary integral and we wanted to fix some variable in it x to be some value, we could fix it with S (x — a) or 
equivalently 6 (/ (x) ) where / (x) is some function which is zero at a; = a. Both delta functions will set a; = a but the 
latter does something else as well, 

Hfi^))-'-!^ (3.1) 

1/ G^)! 

The difference between the two choices of delta function is a factor jy7^^> thus there will be some ambiguity in our 
choice of constraint which we will have to sort out. 

Before plunging into the full field theory calculation with it's infinite number of degrees of freedom (functional 
integrals) it is useful to consider the problem (and solution) in the context of something with a finite number of degrees 
of freedom (function integrals). 
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= y dV e-*^('-) (3.2) 



Consider the above integral, it is invariant under the gauge transformation 6 ^ 6' = 6 + a. Now we would Uke to 
choose a gauge so that we only integrate over gauge inequivalent situations . This means we have to restrict the integral 
so that we only integrate e"'^'-''^ along a line starting at the origin and moving out in the xy plane. Say the function 
= (t>{r) defines this line. (r) should be a many to one function in r and 6 i.e. for any given value of r there is only 
one value of 9 that solves 9 = (j){r). Also, obviously 9 and <f) (r) are defined in the range ^ 27r . In gauge theory this 
corresponds to the requirement that a single physical configuration corresponds to a single gauge field configuration. 
We shall now perform the trick which is central to Faddeev Popov gauge fixing. Clearly we do not change the value of 
the integral if we were to insert a 1: 

1 = j d^r xlx e''^(''\ (3.3) 

Note that we can write 1 = J d(j) (r) 6 {9 — (p (r)) where here one must regard the value of the function at some r, 
(f) (r), as a variable. Regard r as an index for the "variable" ^(r). This means we can write, 

I = /dV [/d(/)(r)(5(6i-(/)(r))] e^^M 

= /rdrd6id(/-(r)(5(6i-(/)(r))e*^M (3.4) 
= J rdr J del) {r) [J d9S (9 ~ <j> (r)) e'^'-''^ . 

Note that though the integrand depends on (r) the integral itself in the square brackets actually does not for no matter 
what value of r we take (j> (r) will always be in the range — > 27r and so the integration over 9 will always kill the 
(j) (r); for every r there is one solution 6* to ^ = (r). This means we can actually shift the integral over (j) (r) behind 
the square bracket, 

I = j rdr j d(j){r) j dOS {9 - (t>{r)) e'^^'''' = j rdr j d9S {6 - ^ {r)) e'^^''^ jd(j){r). (3.5) 
Now for every value of r the integral / d(f) (r) is exactly the same / dcj) (r) s 2tt which gives us, 

1 = 2^ jrdr j d9S {9 - cf) (r)) e*^^''' = / ^^^^ ~ ^'^'^''^ = 27r/0=^(^). (3.6) 

So we have managed to rewrite our integral as a product of a group volume factor (27r) and a gauge fixed version of 
our original integral which will only integrate over gauge/physically inequivalent configurations. 

We shall now perform this example again but in a slightly more convoluted way so as to make better contact with 
the gauge fixing formalism of Faddeev and Popov. The complication is small, we ask the question what happens if our 
gauge fixing condition is in an equivalent but more complicated form e.g. f{r, 6) =0 where the gauge fixing condition 
still has to obey the rules laid out before - one gauge configuration (one 9) for one physical configuration (one r). So 
once again we insert a 1, this time in the form, 1 = J d(l){r)5 (/ (r, 9)), 

I = / dV X 1 X e''^('') 

= /dV[/d<A(r)5(/(r,^))]e'^M 

= /rdrd6»d(/'(r)^(/(r,6»))e'^('-) 

= / rdr / d(A (r) [/ d95 (/ (r, 9)) e'^^'')] 

The next thing to do is to make a change of variables. We know that the constraint / (r, 9) = Q can be written in 
the form/has solutions 9 = (f){r). So the integral in square brackets really has no <p (r) dependence just like before 
when we had instead 5 [9 — (j){r)). Hence we can again shift the integral over J dcj) (r) to the other side of the square 
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brackets where it becomes 27r like in the last example. 

/ = 27r y rdrdeS (/ (r, 6»)) e"^^'') (3.8) 
We must remember to be careful with our delta functions on account of their functional form of their arguments, 

^ 7 = 27r/rdrdi9(5(/(r,6'))e'^M 

= 2^ / rdrd/ (r, 0) (%^) (r,0))e*^M " ^^'^^^ 

Like before we turned our integral into a group volume factor times a gauge fixed version of our original integral. The 
factor ^^^^f^^ now appearing in the integrand typically appears as a Jacobian (determinant) in systems with more 
degrees of freedom, in gauge theory it is known as the Faddeev Popov determinant. In what follows we will see that 
the Faddeev Popov determinant arising in the quantization of gauge theories can be recast in the form of ghost fields. 

So let's see what we can do with QED. 

Z = J DtpDi^DA^" e'^'^ " ''- (3.11) 

Now introduce a gauge fixing condition in the form of a "1". The gauge fixing condition takes the form of a func- 
tional delta function. Such a delta functional fixes the gauge differently at each point in space-time hence we have a 
superscript oo - we have a condition at each point in space-time: 

A (A'*) = j Dx6°- [/ (A'* + di^x)] ■ (3.12) 
Note that our object A is in fact gauge invariant itself: 

A{A^ + ^^'4>) = !Dx5^[f{A^^ + ^^'x + ^^'c|))] 
= /£>x^°°[/(^'' + aM<^ + x))] ■ 

now shift the integration variables x {x) (remember for functional x is like the index, the indexed variables are x {x)) 
at each point by a constant amount given by (j) {x) : 

= n./a;/x(-)'^(^(^) + '^(^)) . (3.14) 
= Y\.xSaiix(:.^)+6(x)d{x{x)+(t>{x)) 

= !D{<P + x) 

Hence we see that 

A (A'' + d'^cj)) = jD{<p + x) 5^ [/ {A^^ + 5^ (x + M 

= ^Dx5'^[f{Ai^ + d^x)] (3.15) 
= A (A") 

is as promised, gauge invariant. We now have a gauge invariant constraint which we can insert into our integral. We 

need to be careful though about inserting this condition, we must make sure we are in fact putting " 1 " into the integrand. 
A (A^) in the general form shown above will non-trivially affect the measure of integration of the path integral. Just 
like in our previous example with r' s and 9' s we want to change the integration variable from x to / (A^ + d^x)^ in 
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doing so we will pick up a Jacobian which we will denote by J. 

= n./d/(x)Det,,,(Mfi)" 
^A(A'') = jDxS^ifiA^' + d^x)] 



(3.16) 



1 



1 



JU,x) 



JDf 6^ [f{A^ + d^^x)] 



(3.17) 



/=o 



X 1 



Hopefully it is clear that the a;, y subscript attached to the Det above means that we are taking the determinant of (...) 
with respect to the x, y "indices". We wish to insert the gauge fixing condition, the "1" in the form, 



l=^(/,x)l/=oAm. 



(3.18) 



The path integral becomes. 



/ DxpD^DA^^ J (/, x)|/=o / Dx 5°^ [f (A^^ + 9^x)] e^^«--' . 



(3.19) 



The next bit of the trick involves changing variables A'^ = + d'^x- Note the action is invariant under 

such a transformation as it amounts to a gauge transformation which would leave it invariant (one can simultaneously 
transform the fermions in the usual way). The other terms are also invariant, we showed above that A [A^^ + d^(p) = 
/ Dx [f {At^ + df'x)] = A (A'') and we will therefore assume that A-^ (A^) = J (/, x)|/=o is also invariant^ 
It then remains to prove that the measure / D-ipDipDA^^ is also invariant, consider for simphcity zero space-time 
dimensions i.e. all the fields exist at a single point xq in space-time. Under the usual transformations A^ ^ A^ + d^X' 
— > e~'-^i/', il) — > e'-^V we have 



/ 



D'il)D'il)DA^ 



I 



I 



D'i\)D'i\)DA>' X Det 



\ 



e^^ 

V 1 y 



/ 



= / D'i\)D'il)DA^. 



(3.20) 



Generahzing to a space consisting of an infinity of points the matrix above becomes an infinite dimensional Jacobian 
which is still diagonal and still equal to one, due simply to the unitary nature of the transformation. Therefore, 



= y Di>Di>DA^^px J{f,x)\f=o 



(3.21) 



So now we have that the integrand does not depend on the gauge x^ this dependence has been absorbed into the 
gauge transformed field Aj^j instead! This means we can bring the integral over x to the front and absorb it in the 



J{f,X)=^et^,y 



Sx{y) 

Now use the functional equivalent of taking the total derivative with respect to a variable, 



J(/,X) = Det:,,j; 



dz 



- Deta: 



(/ 



dz 



sf (Am (^)) 



This object depends on Aj^j but is independent of x- 
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normalization of the path integral as an (infinite) group volume factor {i.e. forget about it) 



iSQ.E.D. 



(3.22) 



leaving the Jacobian (aka Faddeev-Popov determinant) and gauge fixing delta function as the remnants of gauge fixing 
in the path integral. In calculating Green's functions the factor J Dx cancels in the numerator and denominator, recall 
Green's functions are of the form, G {xi, a;„) = — - — ^^^'^^ 



Z[0] SJ{xi)--- 5J{xr,) ■ 

Now we further simplify our situation by choosing our gauge fixing condition such that it is of the form. 
The original path integral has had a "1" inserted into it and after some rewriting now looks like. 



(3.23) 



Z=(^j Dx^ X (^j D^pDx^DA^^^ J{f)\f=o W + Kx])] e^'^^"") • (3-24) 

We should simpUfy this by noting that A'* is a variable of integration, a dummy variable, we should simply rename 

[x] 



Dx]x[ Di,Di,DA^ J{f)\f=o coix) + f{A'^) 



piSQ.B.D. 



(3.25) 



This relation is true independent of what lu (x) is (we didn't specify an u> (a;)), so it will also be true for a linear 
combination of different u {x) 's - with some proper normalization! We will choose a linear combination of all possible 
LO (a;)'s. 



N (C) J Duexp -i J d 



4 W (x) 
X 



2C 



1 



(3.26) 



where N {Q is the normalization we just mentioned. Therefore in our linear combination of Z's each one is weighted 

2C 



by a factor N {(,) exp — « J d a; '^^c '•^^ normalization N {Q) is making sure that when we add them all up we 



get just Z again. Hence, 

Z = N{Ox{jDx)x 

J Doj&xp -ijd'^x^^ (^J D^jJD^DA^' J(/)|^^o uj{x)+f{Af') 
= N{C)x{jDx)x 



piSQ. 



exp — i J d'^x 



4 i^jxy 

2C 



/ DuJDi}Di)DA^' 
= 7V(C)x(/Z)x)x 

QDipDipDA" J(/)|^^o e*^Q-^°exp -i/d^'x^) 



^(/)l/=o u;{x)+f{A>^) 



oiSQ.E 



(3.27) 



In the last step we have integrated over oj (x) which has removed the delta functional and turned exp — i J d*a; "^^^ 
into exp — i J d^a; It is worth taking a breath here and making sure you are happy with the above which may 

look like smoke and mirrors, it is not. Remember what we said above about the Physics (Green's functions) being 
independent of the normalization of the path integral, consequently we will disregard the normaUzation, instead we 
will work with Z = ij^y 

To simplify further we need to actually specify a form for / (A^). The form that is commonly used is the so-called 
covariant gauge (it is clearly Lorentz invariant). 



f{A^) = d'^A^. 



(3.28) 



The gauge condition above is such that if F^^" is something then F'*'' is something else, where F^^'^ is the same as F^^" 
but with the replacement, A^ ^ A^ + d^cf). So what does this make the Faddeev-Popov determinant? 
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•^(/>x)l/=o = Det,;,y 







Sx(y)) 








Sx(y) 



(3.30) 



'■x,y 

= Det,,, (9^9/^(5 

In the Faddeev-Popov gauge fixing formalism the determinant is rewritten as a functional integral over an exponential 
of ghost fields with the operator above sandwiched between them. 

A Grassmann variable / number is an anticommuting number. For a set of Grassmann variables 9i this means, 
{6i, 6j} = 0, this trivially means that the square (and therefore also higher powers of a Grassmann variable) equals 
zero as anti-commutation requires 9i6i + 6i6i = 0. Consequently Taylor expansions of functions of Grassmann 
variables temiinate after a few terms. These peculiar "numbers" also generate ambiguity in the definition of integration 
over them. In fact the most natural definition of integration for Grassmann variables turns out to be such that it is 
the same as differentiation. We attempt to bolster these statements in the last section (5.3), for a dedicated discussion 
of Grassmann numbers we refer the reader to [17], for now we will merely state the definitions. Taylor expanding a 
function / {9i) of a Grassmann variable 6i we have, 

f{9i)=a + bei. (3.31) 

Integrations over a and b9i are defined as (a and b are "normal" bosonic numbers). 



(3.32) 



/ de, a = & J- (a) =0 
jd9ib9i = b & 4:{b9i) = b ■ 

Jd9if{9i) = b 
also 

(3 33) 

Idic9,)a = -^(a) = 

Id{c9i)b9i = 3(4y(i6ceO = -cb 

An important point to take away from the last two equations above is that for Grassmann variables the integration 
measure changes in the opposite way that it does for normal variables i.e. 

J d {c9i) f {9i) = ljd9if {9i) . (3.34) 

Also, for complex 9i we simply treat the two components (real and imaginary) as two independent Grassmann variables. 
Finally, as you may well variables and therefore also the integrations, anticommute. 
Now consider the following integral, 

]^ / dwzdrjz exp i^^'EyTja^T^yLOy (3.35) 

where x and y are discrete indices of complex Grassmann variables r]x and Wy. is understood to be the complex 

conjugate of lo more commonly denoted in the literature by a bar over lu: rj ^ uj) . T^y is some Hermitian matrix of 
("normal" / "bosonic") numbers. We would like to evaluate this integral. We don't know what to do with it as it stands 
assuming T^y is not diagonal. First we want to transform the integrand and variables of integration such that the F in 
the exponent is diagonal. To do this we introduce a unitary transformation matrix U as follows, 

^x'^V Vx^xyi^y = '^x'^y'^a'^b^c'^d flxUljJab^bcUljJdy'^y (3.36) 
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We choose U such that T,bT,cUab^bcUl^ = fad is diagonal i.e. fad = ^aa^ad^- Note that the we show the summation 
over X and y expUcitly with E, a repeated index does not imply a sum over it. Now we redefine our fields according to 
the unitary transformation such that 

'^x'^y Vx^xyi^y = '^a'^d Va^aa^^a^ad = '^a Va^aa^^a (3.37) 

i.e.,s 

Rewriting the integral 3.34 in terms of the transformed variables we have that, 

exp iY^x^yVx^ xy^y — 

Yl d [EpUlgQp) d i^gfjgUp^) exp iE„^„f„„w„. (3.39) 

We need to get the integration measure in a more friendly form. As we said before, the Grassmann integration measure 
transforms as the inverse of the Jacobian of the transformation instead of just the Jacobian as is the case for regular 
numbers. For a slightly better proof of this consider the following integral over complex Grassmann variables, 

j AOn- j d6»2 j AOi 6162. ..On = 1 (3.40) 

as J d6i 6i = 1. Now we make a change of variables 9i T,jUij9j. The change of variables cannot affect the value 
of the integral. How does the integrand change? As the ^j's are Grassmann we can write. 



'iiOi^.-.Oij^ — —Y^j^Y,j^...lljj^e^'^^^'"^'^6j^9j^...6jff (3.41) 



i.e. we sum over all A^! permutations of the A^ variables remembering to divide by A^! at the end and keeping track of 
the minus signs with e-'^-'^ - -'" . Therefore we can write the transformation of the integrand as. 



?i^2-"^iv 6x62. ..On = -^EjjSi2...Sj„Sj^Sj2...Sj„e-'^-'^'"^^^ji^i2...^,^ Ui^jiUi2j2---Uiffjfj . (3.42) 



We can use the anticommuting property again to simplify further and write that for a given combination Oj^ . . .6j^ , 



e'^'^-'''9i92...9N (3.43) 



note that here we have no implied summation, there is no sum taking place above and there are no indices on the thetas 
on the right! 

9\92...0n — > 9i92...0n 

= (^E,,E,,...I]i,E,-,S,,...S,,e^i^-^«e^i*--" Ui.j.Ui^j^-Ui.j,) 9,92...9n 



Now the (definition of the) determinant ofanNxN matrix M is. 



m 

(feel free to check this) so. 



Det(Mi,) = -i-E„,...E„„E^,...E;3«e""-""e'^"-'^"M„,^,...M„„^^ (3.45) 



9i92...0n kh-ON = (Det {U)) 0i92...0n. (3.46) 
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Consequently if we need, 

d0i 9i92---0n = J d9N:.J d92 j d9i 9i92...9n = 1 (3.47) 

then we must have that, 

J ddN... j d92 j dh = :^-L_y d9N... j d92 j d9i. (3.48) 

So invariance of the value of a Grassmann integral under a simple change of variable means that the Grassmann 
integration measure transforms as the inverse of the Jacobian of the transformation as opposed to just the Jacobian for 
regular numbers. 

In our case this means that, 

n / dc^.n / - Bet (m) Bet (m il / ^-^11 / d% = R / d-^II / d% (3-49) 

X y ^ ^ ^ ^ X y X y 

as WU ~ I. We can manipulate this equation above by moving the ?7's w's and their corresponding 77's and oj's in the 
same way on both sides of the equation without worrying about minus signs. This is because every time we reorder 
things on the left side we may generate a sign but the same reordering on the right wiU naturally generate the same sign 
there so we can forget about signs and say, 

n / du;,dv. = p^t(^t|Bet(t/) H / ^^'^^i^ =11/ ^^^'^^^^ (3.50) 
The result is that the entire integral is invariant under the transformation which diagonalizes T^y, 

Yl / dw^d??^ expiT.a:T,yr]a:Ta:y(Jy = 11/ d'^^^^^ exp iSa^afaaWa (3.51) 
z z •' 

Taylor expanding the integrand on the right gives, 

j dUzdr]^ expiT,x^yT]xrxyOJy = H / '^^^^^ (^^ + ^a??afaaWa + ^ (^^aVai'aal^a^ (pbVbi' bbl^b^ + 

(3.52) 

For a term in the expansion to survive the integration it must have one copy of each of fj^ and Cjy for all elements in the 
products of integrals. Take for example J duzdfjz = J dujzidfizi J dCjz2dfiz2 i-e. there are only two values of each 
index. The expansion terminates after the second term (or more generally for N values of the indices, the A'^th term), 

exp iT,afiataai^a = 1 + ^a^of oa^o + ^ (^aVai'aal^a^ {^bVbi' bbO^b^ (3.53) 

as higher terms will have to involve the square of at least one of the Grassmann variables. The first term (1) vanishes as 
J d^ 1 = for any Grassmann variable by definition and we are integrating over four different Grassmann variables. 
The second term also vanishes for basically the same reason. This term is simply, 

Cox2 (3.54) 

but the integral involves integrating each of these over four different Grassmann variables so they both vanish. Only 
the last term makes a contribution to the integral. It contains terms with one r] and w for each r] and a; integration. The 
last term is 

^X\^X\X\^X\Vx2^X2X2^X2 ~^ 'nX2^ X2X2^X2V^X\^ X\X\^X\ ~ {jix\^X\fix2^X2 ~^ 'nX2^X2Vxi^Xi)^ X2X2^ X\X\ ■ (3.55) 
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where we have fjxii'xixi'j^xiVxii'xixii^xi = as they contain squares of Grassmann variables i.e. the last term is a 
sum of the cross terms between the products of sums. If we had N values of the index it would be the A^th term in the 
expansion that would survive, for the same reasons. We are free to split these terms into pairs of adjacent Grassmann 
variables and freely move these pairs around without worrying about picking up minus signs, such pairs commute with 
each other, 

[^1^2) ^3^4] = O162O364, — 63646182 = 61626364 — 61626364 = (3.56) 

which basically means all of the aforementioned cross terms in the integrand are the same. How many such cross terms 
are there? Easy, we need a different rjTijo from each sum, the A/'th term is a product of N such sums and each sum has 
A'' different jjEw's in it. So I can get the non-vanishing term N\ times corresponding to the number of permutations of 
the N pairs. 



= Dz/dwzrfjy^ M'72ifziziW2i...?7z„f2„z„Wz„ (3.57) 

= (riz f (riz / dujzdf]^ ^zjv'^^iv-»?2i'^2i) 

As r is r diagonahzed this means that the first product above is the product of the eigenvalues of F which is the 
determinant of F: 

Det(fzz) = n.f.. 

= Det(C/rC/t) ^333^ 

= Det(C/)Det(F)Det(?7t) 

= Det(?7t?7)Det(F) 

asDet(Z7t?7) =Det(/) = 1. 

l\^Jdu}:,dr]zexpiT,x'Eyr]xTxyUy = (Det (F)) (H^ / dwzd??^ Tyz^w^^.-.Ty^jW^J 

= (Det(F)) (/dwzid7?z,.../da;z„d7?2„ 7?zj^Wzjv-"%i'^zi) (3-59) 
= (Det(F))xl 

Generalizing to the continuum limit we have, 

j DcoDrj exp i j ^^^^ j d*2/ V (x) F (x, y) u (y) = Det^,^ (F {x, y)) (3.60) 

and finally we have a way of including the Faddeev-Popov determinant in our path integral in terms of "Ghosts" r] and 
LO. In the case of our covariant gauge we have. 



'DeXx,y {d^di^S {x -y)) = / DuiDri exp i / d^'x / d!^y rj (x) d^d^'S {x-y)uj (y) 
= J DloDt] exp i J d'^x rj (x) d^d^u: {x) 



Finally the gauge fixed path integral is, 

Z = (0|0) = j DA^'D^D■^l)DuJDr] exp i j d*x j-^F^^i^^"^ + ^ {d^A^'f + rjO^d^u + P-m)il>^ 

(3.62) 

with a kinetic term rjd^d^uj appearing in the original QED action due to gauge fixing. The ghosts are not coupled to 
any of the physical fields and so they don't appear in perturbation theory. In non-Abelian gauge theories fixing the 
gauge results in ghosts as in QED but the resulting action will have a coupling of the ghosts to the physical fields. This 
couphng results in the ghosts appearing in Feynman diagrams, in perturbation theory. The ghosts are so-called because 
they are unphysical for a number of reasons, the most obvious one being that they are Grassmann fields yet they are 
also (complex) scalar fields (pseudoscalar particles) they have no spin! This unphysical spin- statistics relation means 
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that ghosts cannot appear as external particles / external legs on Feynman diagrams, the perturbation expansion is such 
that (in non-Abelian theories) the ghosts only appear in closed loops. 

Finally, as ghosts are an artifact of gauge fixing they manifest themselves in the Lagrangian and hence in perturba- 
tion theory in different ways according to the choice of gauge, in certain gauges the ghosts may not even materiahze 
in perturbation theory. However gauge fixing affects the photon propagator (see next section) and in most instances 
it turns out that calculations in perturbation theory are made much easier by choosing a gauge which simpUfies the 
propagator of some physical particle at the price of including ghosts. 
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3.2 Feynman Rules in QED. 

We now have a fully gauge fixed path integral for QED ghosts e/c...(put in fermion fields now). 



Z = (0|0) = y" DA^'D'^D'^Dt]Duj exp i J d^x i^-^Fi^^F"" + ^ (S^A")^ + rjd^'di.w + ^ (i /) - m) v| • 

(3.63) 

Firstly we will try and derive the photon propagator, we therefore want to calculate the following (we are normalizing 
everything by 



(0|A''(a;i)A''(x2)|0) _ J DA''...Du A'' {xi)A-' (x2) exp iSQ.B.D 

(o|o> - Z 



5M.f)5M^2) 1^ / DA^'-Di^ exp i {Sqed -ijd^x J« {x) A« {x)) 



57^(xi)(5J„(x2) 



..=0- (3.64) 



We have defined Z [J'^] = J DA^ ...Dlo exp i {Sqed — if d^x J„ (x) A'^ (a;)) with J" (x) a source term. The source 
term allows us to rewrite our definition of the two point function as above, it is of no real physical significance. Though 
J'" (x) is introduced ad hoc into the Lagrangian it is not affecting the Physics as can be seen by the equivalence of 
everything in the working above. The next thing we do is make a perturbative expansion in e. The interaction part of 
the generating functional is expanded in powers of the coupling constant. 

exp j d'^xe^) /A^ = l + e j d'^x M - O (e^) (3.65) 

It is also necessary to do some rewriting of the free part of the Lagrangian. Also up until now we have not had to worry 
about the (space-time) variable that the fields depend on, we have safely been able to assume everything depends on x 
but things get a bit more complicated here. Consider the bit for the photon (the same procedure works separately for 
the fermions). 



jd^x-\ (d,A (x)^ - d.A (x) J {d'^A'^ (x) - d-'A^ {x)) - ^ {d,A>^ {x)f - J. {x) A [xf ■ (3 



66) 



For reasons that will become apparent the next thing we want to do is to put this in a form J d^x A^Q^^^A^. Forgetting 
about the source term we have 



/ 



d^x - i {2d^,A, {x) {x) - 2d,A^ (x) A' {x) + -^d^^A^ (x) d.A" (x) ) . (3.67) 



Integrating by parts gives 

= / d^x - i (-A, (x) d^d^A'^ (x) + A^ (x) d'^d^A, (x) - \A>^ (x) d^d^A'^ (x)) 

\ ^ ' . (j.oo) 

= / d^x - i (Aa {x) {-g'^Pd^^d^ + (1 - C"') d^d") Ap (x)) 

We now need to sort out the source term which we so conveniently forgot about, this is done by shifting our field 
A^ (x) everywhere such that 

Ai^ [x) ^ A'" (x) = Ai^ (x) + j d^y A>^^ (x - y) J. (y) 
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where Q**^ (a;) = iSj^S"^ {x) i.e. Qi^" J d% {x - y) J" {y) = Ji^ {x), this gives, 

-lA^Q^^-A, - J,A- ^ -i (a; (x) - J d^y (x - y) J- (y)) Q^-^ (A', (x) - /dS (a; - z) (z)) 

-J, (x) A'- (x) + J. (x) / d^y A-^ {x - y) J, {y) 
= - \A!^ {x) Qi^-A', [x) + i a; {x) Q>^- f d^z A^x {x - z) [z) 
+ i (/ d^y A^« {x - y) J« (y)) Q'^M', {x) 

- \ (/ dS A^« {x - 2/) J'' (y)) Q^- (/ d^y A,^ (a; - z) r {z)) - J« {x) A!^^ [x) 

+ J, (x) / dSy A-'' (,x - y) J„ (y) 
= - i a; (x) g^M^, (x) + ^A'^ (x) (x) + i (/ dSy A^, (x - y) J« (y)) Q'^M', (x) 

-i (/ d^y A^, (x - y) J« (y)) (x) - (x) A'^^ (x) 

+J«(x)/d^y A«- (x-y)J,(y) 
= - \A'^ (x) Q^M', (x) - i a; (x) J'' (x) + i (/ d^y A^« (x - y) J« (y)) Q"-- A', (x) 

+ ij«(x)/d^y A^-- (x-y)J,(y) 

(3.69) 

A long winded and tedious calculation^ involving the use of the product rule and setting surface terms to zero means 
that the 3rd term can be rewritten, 

i (^j d% A^« (x - y) r (y)) Q^^A!, (x) = \ {q^^ j d^y A^^ (x - y) j' (y)) A!, (x) (3.70) 
By definition of the function A^^ (x — y) from above, this equals 

\r{x)A'^{x). (3.71) 
So to first order in e the generating functional can be written, 

Z[J]k j DA!^'D%PD%PDt^Duj (^ + e j d%i> sxpiSpreeQED (3.72) 

where 

SpreeQED = / d^x / d^y i J„ (x) A^^ (x - y) (y) 

/ d*x - i a; (x) Q^M', (x) + 77 (x) d'^df.ij (x) +i;{x){i jd-m)iP (x) 

i.e. QED without the interaction term and the source and field terms rewritten. 

We now arrive at the conmion sense result that to get the photon propagator, the photon two point function, we will 
be taking the zeroth order of perturbation theory i.e. ignore the interaction terms. If we take higher orders the higher 
order terms in the interacting part of the Lagrangian will be playing a role in our correlator i.e. we would be having 



i (/ (x - y) J- (y)) Q^^-K i^) 

= i (/ d*y A^s (x - y) {y)) {-g^^-d-d^ + (l - C"!) 9" 9") K W 



la- ((/ A^s {x - y) (y)) i-a'^-d^Al (x))) - {^d^ f d^ A^s (x - y) (y)) {-g'^^d^A^ (x)) 
ia" ((/ d-'y A^s (x - y) J« (y)) ((l - C'^) d^A', (x))) - i ((9^ / d^y A^g (x - y) {y)) ((l - C'^) K {x))) 
\d^d- {{J d-^y A^s {x-y)j'{y)) {-g'^^'K (^))) " ifl'' ( (9« / d^y A^^ {x - y) {y)) {- g^^- A',Ax))) 
-5„ / dS A^s (x - y) (y)) {-g'^^K {x))) + {Id^d- f d^y A^s [x - y) (y)) {-g^^-A'„ (x)) 

+ ^df^d- ((/ d^y A^s (x - y) J* (y)) ((l - C'^) K {x))) - ((9" / d^y A^s (x - y) J' [y)) ((l - K (^))) 
\d- ((aA' / dS A^s (x - y) J' (y)) ((1 - C-^) K (x))) + \ {{8^8- J dS (x - y) (y)) ((l - C'^) K (x))) 



Assuming that the field A^^ (x) and its first derivatives vanish as a; ^ oo we can remove six of the above terms on the grounds that they can be 
revratten as surface integrals. 

= k i-g'^'d^a'^ f d^y A^s (x - y) (y)) K {x)+\{{l- C'l) di^d" J d^y A^e {x - y) (y)) A', (x) 



2 



(Q"" / d^y A^s (x - y) (y)) A', (x) 
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terms with e / d^yip Ai^ bits in them in our correlator. So to get the photon two point function {OlA^^ {xi) A'^ {x2)\0) 
to zeroth order in e (normahzed by Z~^) we want to functionally differentiate, 

Z[J]= j DA^Di^DijDrjDu (1 + O (e) + ...) exp iSpree q.e.d. (3.74) 

with respect to J^{x\) and (2:2), 

(0|A^'(a;i)^"(a;2)|0) _ S'^ \n Z[.J] 

(0|0> - S.J^{Xl)SJ^{X2) 

= / DA^D{^D^PDtjDlo A'^" (xi - X2) exp iSpree Q.E.D. + O {e) + ... (3.75) 
= ^((0|AA"'(xi-X2)|0)+O(e) + ...) 

So the propagator (to zeroth order in e) A^" {xi — 0:2) is essentially just the inverse of the differential operator Q^^, 

Q^'-A,, (x) = {-g^'d-'d^ + (1 - C"') d^d"-) A,, (x) = tS^iS^ {x) (3.76) 

Fourier transforming we have, 

{-g^^-^dW^+{l-C-')d^^d'')Jd^kA,^{k)e'''-- = iS^S^ (x) 
/ d^'fc A^« (fc) (^^-^fc^ - (1 - C"^) fc^r) e**^-^ = iSaS^ix)' 

So now the question is how do we solve A^.^ (fc) {g^"'k'^ - (l - C~^) k^k") = iSj^ for A^« (A;)? To start with let's 
rewrite Q^P; 

^ - I 5 + fe. J c- fe^ (3.78) 

_ i.2 pAp fe2 pAp 

— ft r-y T ■ 

Where we have defined the projection operators, P^^ = g^^ — ^-j^ and P^'^ = -^^l-. As one would expect PtPt = 
Pt, 

~ y /s2 fe2 (-3 79-) 

- 9 p- 

_ pA^i 

— I- J, . 

The same is true of Pl, PlPl = Pl- Also PlPt = = PtPl, 

P^Pn P***" — fc^fc" „ („Kp _ k'^k'^ \ 

gpK-^T — fe2 i/pft \g 1.2 ) 

= . (3.80) 

= 

This last identity shows us that Pt and Pl project out orthogonal subspaces, i.e. the Hilbert space of states is now 
split into two by Pt and Pl, all the states in one half being orthogonal to all the states in the other. 

Returning to the derivation of the Feynman rules we recall that we are essentially after the form of the object Q~^. 
Given our study of Pt and Pl and their relation to the space of states we will assume that can be written as a 
linear combination of these two operators: 

Q-i PP- = ^pw + (3 81) 

Naturally one QQ~^ will give the identity or equivalently Pt + Pl- 



Q^'gp.Q-'^" = -k^ [P^' + ^^P^') gpAAPT^" + BP^^) 



pA^ I pA^ 



(3.82) 
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(^AP^" + f P^^) = P^" + P^ 




(3.83) 



Given the form of the photon propagator one might wonder why then we don't see ('s in our Feynman rules for QED. 
The answer is simply that the photon has no longitudinal component. This is not to say that if we were to put our 
photons inside some pathological loop diagram we would not see the ('s, we would but they all cancel in the final 
matrix element. The proof that all the C's cancel is not obvious and we come back to it in the next lecture when 
discussing the QED Ward identities. 

The same recipe can be used to give the following propagators for the ghost fields and the fermion fields... 



CO 



► 



l/—m+ie 



Note that the e's in the above Feynman propagators are put in by hand to regulate the poles in the Feynman propagators 
(branch cuts) where it is understood that in the calculation of the amplitudes the limit e ^ is taken at the end. 
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3.3 BRS Symmetry. 

Continuing from lecture 7 - the generating functional for QED with ghosts ere... (put in fermion fields now). 



Z = (0|0) = ^ DA^'D'^D'^Dt]Duj exp i J d'^x l^-^Fi^^F"" + ^ (S^A")^ + rjd^'di.w + ^ (i /) - m) v| • 

(3.84) 

This thing looks gauge invariant except for the d^Af^ gauge fixing term. It turns out that despite the apparent breaking 
of gauge invariance above there is in fact some residual symmetry present, we shall come to this point shortly. This 
symmetry was the result of work of Becchi, Rouet, Stora and (independently) Tyutin, it is of crucial importance in the 
quantization of gauge (and other more general) field theories. The BRS(T) symmetry is based on the notion that we can 
have a further symmetry of the action if the ghost fields transform non-trivially. BRS is a generalized gauge invariance. 
In general the parameters of gauge transformations are arbitrary position dependent functions. In the BRS arrangement 
this arbitrary function is considered as being a product of two Grassmann quantities, one a constant the other a function 
of space-time. Consequently anything that is gauge invariant is BRS invariant - i.e. the original action, prior to gauge 
fixing, is gauge and therefore BRS invariant. It has been the nature of these notes to exhibit such obvious results by 
brute force too. 

Suppose that under these new BRS symmetry transformations we have: 

A^^A'^^A^ + ed^Lu (3.85) 

where this e is an infinitesimal Grassmann parameter and lo is our ghost field from before (also Grassmann). Under 
this transformation the field strength tensor is in fact invariant: 

p . fpi f) A' — Pi A' 

= d/^A^ - d^Aij, + edij,d^u} - ed^d/^uj 
— 9fj,Ai, — di,A^ 

What happens to the gauge fixing term under BRS? 



-^^{d^A^^f-^{d^A^){ed,d-'u;) 



Well it changes, so how do we fix this? We must specify a way for ghost field 77 to transform so as to cancel the 
— ^ (edud^uj) coming from the transformation of the gauge fixing term in the action. We do this by making 77 

transform as 77 ^ + ^ (i9^A^) e then the ghost term in the action {rjdvd^u) turns into itself plus a part which cancels 
the aforementioned term from transforming the gauge fixing term. In addition to this we make omega transform into 
itself w — > w and the fermion fields transform as — > e^^'^'^tjj, -tjj — > e~^^'^'^'ip. It is worth Taylor expanding the 
exponential in the last pair of transformations. 

,1 s (3.88) 

= (1 + ieeu) + 5| {ieeu}) (ieeu)) + ...) tp 

Note that e and ui are Grassmann numbers so = etc which means that the above Taylor series 3.88 terminates after 
the first two terms! 

^ ^ -)/,' = (1 + ieeu) tp (3.89) 

Likewise we get for -ip, 

^ ^ tjj' = {l-ieeu;)xp. (3.90) 
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So for the fermion part of the Lagrangian we have, 

'ip {i /J — m) ip i/j' (i /D' — m) tp' 

= (1 — ieeto) {i /) + e ^ + ee (jduj) — m) (1 + ieeco) tp 

= (1 — ieeoj) "ip [{i jd + e /i)ip + ee [jduj) ip — mip — ee jd {ojip) + ie^eoj /itp — imeeojip) 
= '4) {{i jd + e /A) ip + ee {jdui) xp — rmp — ee jd (cj^j) + ie^ecj /iip — imeecutp) 

— ieeui'ip [{i /) + e /A)^j + e^e [ jdu) i]j — milj — ee jd (ujip) + ie^eu) /Atj) — imeeLOip) 

= ^ {{i jD — m)ip + ee {J)uo) ip — ee jd (i^'tp) + ie^e /Aujip — imeeojip) — ieeu'tp {i — m)ip 
= ip{i fi — to) Ip + eipe {fiijj) Ip — eipe f) {ojip) + ie'^ipe /Acoip — imeipe {ojip) 

— ieeu-ip {i fi — rnjip 

(3.91) 

Plug in: —eipe jd (cotp) = —eipe{jdu}) tp — etpeu) {/dtp) ■ 

= -ip {i /D — m) Ip — e-tpeui {/dtp) + ie'^-ipe /A {ujtp) — imetpe {totp) — ieeco-ip {i /D — m)ip 

= "ip {i /D — m) tp — etpeu {/dtp) + ie^tpe /A {wtp) — imetpe {u)tp) — ieewtp {i jd + e /A — m)tp 

= tp{i /D — m)tp — etpeu {/dtp) + ie^tpe /A {utp) — imetpe {utp) + eeutp {/dtp) — ie^eutp /Atp + imeeutptp 

(3.92) 

Now if we shuffle the variables around in the 2nd 3rd and 4th terms around taking into account the anticommuting 
nature of the Grassmann variables tp,tp,u and e we have, 



—etpeu {/dtp) = 



+ie'^tpe /A {utp) = 



—imetpe {utp) = 



+eetpu {/dtp) 
—eeutp {/dtp) 

+ie^tpeu /Atp 
—ie^etpu /Atp 
+ie^eutp /Atp 

+imeetputp 
—imeeutptp 



(3.93) 



(3.94) 



(3.95) 



(3.96) 



=> tp{i fi — m)tp tp' {i fi' — m) tp' 

= tp{i — m)tp 

We have a new symmetry of our gauge fixed QED action. Is BRS also a symmetry of the path integral measure? In 
other words is the Jacobian between the variables Ai^ {x) ,tp {x) ,tp {x) ,7] {x) ,u {x) and their transformed counterparts 
equal to one? Recall the BRS transformations. 



{x) - 


^ A'^'{x) 


= A'' {x) 


+ 


ed'^u {x) 


tp{x) - 


tp'{x) 


= tp{x) 




ieeutp {x) 


tp{x) - 


tp'{x) 


= tp{x) 


+ 


ieeutp {x) 


r]{x) - 


T]' (x) 


= v{x) 


+ 


^{d^A^^{x))e 


u {x) - 


u' {x) 


= u {x) 







(3.97) 
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Imagine space-time consists of a single point (a;o), the Jacobian would be of the form: 

/ dA'^(xa) dtpixp) d^(xo) dri{xQ 



Det 



Det 



<?Af (no) 
aA"(rto) 
di:(xo) 
dAi^(xo) 
dip{xo) 

dr}(xo) 

du>{xo) 
1 



dA_^{xo) 
dtpixp) 
dtpixo) 
djiixo) 
dtjj(xo) 
dipjxo) 
dnixo) 
dijiixo) 
du){xo) 





dA>-(x„) 
d^(xo) 
dtpixo) 
dip(xo) 
dt/j(xo) 
d-ipjxo) 
dri{xo) 
di>{xo) 
dui{xo) 



dAi'(x(,) 
dr^jxp) 
dip(xo) 
drijxo) 
diji(xo) 
dyixp) 
dri(xo) 
dTi{xo) 
dui{xo) 



9oj{xa) 
dAi-ixo) 
<i^(xu) 
dip{xo) 
9u}{xa) 
dip(xo) 
duijxa) 
dri{xo) 
duj{xo) 
9w(xo) 







1 — ieeoj (xq) 

1 + ieeu) {xq) 



\ O (e) ieeip (xq) iee^ (xq) 

Due to the nature of e the off-diagonal elements do not contribute so 

J 



0{e) 


\ 














1 








1 ) 



(3.98) 



= 1 X (1 - ieeuj {xq)) x (1 -|- ieeco {xq)) x 1 x 1 
= 1 — ieeoj (xo) + ieeco (xo) + e^eco {xq) eoj (xo) ■ 
= 1 



(3.99) 



This is exactly what we want. In reality the Jacobian takes the form of the determinant of an infinite dimensional 
matrix (because there are an infinite number of space-time points in real life) but hopefully it is still clear that it will 
still essentially be like the matrix above repeated down the diagonal an infinite number of times (one for every space- 
time point). The resulting determinant will appear just like above with only terms linear in the uj (because all other 
terms higher order in iv will have at least an attached to them and the linear terms cancel each other out at every 
space-time point just as they did above. 

The partition function is totally invariant under BRS, so we can start generating Ward identities in the usual fashion. 
Consider the general case where we have some time ordered product of fields given by the function H [A^, -tjj, tp, r], ui) : 

{0\T\H {Ai',^,^P,Tj,Lo))0 = J DAt'Dij}DijDr]Du H {At',iP,ij,r],uj) 

exp I / d^x (i - m) V - \F^.F^''' + ^ {d^At' f + 7?a^5^w} ' 

(3.100) 

In not specifying the form of H we are essentially studying the most general type of Green's function. Now we do the 
usual change of variables trick on the above (BRS transformation): 



A^'{x) 
tp{x) 

il;{x) 
r]{x) 
u) (x) 



A'^" (x) 
i,'{x) 

T]' (x) 

uj' {x) 



Ai^ [x) 
tp{x) 

il^ix) 
■q{x) 
u> (x) 



+ 
+ 



ed^LJ {x) 
ieeoj-il) {x) 

ieeujij) [x) 
\{d^A>^ {x))e 



(3.101) 



We know that DA^'D^pDipDrjDu and exp i J d'^x ^t/j {i P - m) tjj - \Ffj,^F>"' + ^ {^^A^'f + rjd^'d^u^ trans- 
form into themselves from our study above but this need not be the case for H, it isn't necessarily BRS invariant. So 
what happens to HI H ^ H' = H + S^H. 



5,H 



M^S.A^ + ^6,v + If <5eVi + If ^.V' + ^S,^ 



^ed'^uj {x) + %6,n - ie^eoji, + ieff ec^V + %Q 



dAi* " 
dH 



dri 
dH ; 



^ed^uj{x) + ^5,r, 



«e||ewV; -I- ie^euji) 



(3.102) 
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Under the (BRS transformation) change of variables we have 



iS 



+ J DAi'DxpD^DriDLO {S^H {A'',^,iP,'q,uj)) e*^ 

Using Dirac notation and substituting in for ScH we have, 

{0\TH{A^,i;,i;,r,,uj)\0) = (O |T F (^'*,^, V,?7,w)| O) 

+ (0\T6,H{A^',^,tP,rj,u;)\0) 



(3.103) 



(3.104) 



i.e. we have the Ward identity 



^ 9ff ^„ , , dH ^ . dH - . dH - 



> = 



(3.105) 



Note that this Ward identity does not depend on e, we have not required the parameter of the BRS transformation to be 
small! The Grassmann nature of e terminates any power expansions after the linear term appears (the same reasoning 
applies to the case of non-Abelian theories which is a simple extension of our formalism above, instead of e we have a 
vector of e's, one per gauge group generator). 
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3.4 BRS in QED. 

Recall the BRS transformations from the last lecture: 



Af" (x) - 


A'f" {x) 


= At" {x) 


+ edf^uj (x) 


ip{x) - 


4>' (x) 




— ieeuj-ip (x) 


tp{x) - 


V' (x) 


= 


-f ieetotp (x) 


r]{x) - 


r/ (x) 


= 


+ ^{d,A^{x))e 


u (x) - 


J {x) 


= U) (x) 





(3.106) 



We also discovered that the BRS symmetry gave us the following general Ward identity for some string of fields H: 



(0 \T6eH\ 0) 



^ dH , ^ dH ^ . dH - . dH - 



) = 



(3.107) 



needless to say these BRS Ward identities can get arbitrarily complicated, so let's consider an easy one first. Consider, 

H = r]{x)A^{y) (3.108) 



5,H = r,{x)edjF^uiy) + ^edi''^A-{x)A^iy) 
= e l^^di^^A'^ (x) A^ (y) - r, [x) djT^u; (y)) 



(3.109) 

Note +r] {x) ed^uj [y) has become — er; {x) (y) on the grounds that e and rj are Grassmann numbers. Differentiation 
with respect to x has been denoted 9^^^ and differentiation with respect to y has been denoted dl^^ etc. So we get our 
first BRS Ward identity, 

/niTia^^^ 4" (x) A.. (7A m = (n\Tn (x^ a'.^L (v) m 

(3.110) 



{0\T^dl-'>A'^{x)A^{y)\0) = {0\Trj (x) d)I" u; {y) \0) 
di''\0\T'^A'^{x)A,{y)\0) = djl'\o\Tri {x) uj {y)\0) 



which relates the divergence of the photon 2-point function {OlT^A" (x) A^ (y) |0) to that of the ghost two point 
function (0 ITr] (x) u) {y) \ 0). Recall that in QED the ghosts do not couple to the other fields so we already know what 
it is! 




fc2 + ie 



9p^ - -p- (1 - 



fc2 + ie 



{PtpX + CPlpx) 



(3.111) 



11 
I 



CO 



fc2 + ie 

Where PrpX and PlpX are the transverse and longitudinal components of the photon two point function. 



(3.112) 
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We shall now attempt to simplify matters by Fourier transforming our Ward identity. For the ghost part we have, 

d'tl'\Q\Tn{x)u^{y)\Q) = F.t\^) 

— r d-'fc ik(x-y) 
J (2^ k-'+ie'^ 

For the time being we feign ignorance of the photon propagator in the hope that being general may afford us some way 
of incorporating modifications to to it by e.g. fermion loops, hence we denote it by which means the photon part 
of the Ward identity transforms as: 

Plugging this into our Ward identity we find. 

There exists a theorem which states that if the Fourier transform of something is zero then that something is also zero! 
Hence, 

Vi)^.(fc) = -^. (3.116) 



Now let's write the 2-point function as follows: 



♦ 



ry = i:}(0) -I- (Ot-jP^^ (O) 



where D^^} is the free photon propagator given above. In our Ward identity we have ^k'^D^^ (fc), substituting in the 
above this gives 



ik-^D^^ {k) = f fc^ { {Pt,. + CPl,.) + {k) W>^d'^^1 (fc)} 

k^-\-ie ' 

Use the definitions of Prpx and Plpx (i.e. in terms of fc's) to simplify the above first term inside {...} 

f^^-feq^iJ ((s/ii' fes^) + C ("P^)) = ^■fe2+jj {{kp. - kp) + C,kp) 



{k'^D^^ik) = -^ + ^r4°;(fc)n^^i)i°J(fc) 



(3.117) 



(3.118) 



Substituting this into the Ward identity on the previous line gives 

k-'Dl^^ (k) n^^Dfj (k) = 0. (3.120) 

This already looks pretty interesting as D^p' (fc) HP^uf^ (k), which was a generic way of writing any and all radiative 
corrections to the photon propagator, is looking heavily constrained if not totally constrained as a result of our simple 
Ward identity. So let's take a closer look. Let the fc" hit the free photon propagator D^^ (k) (k" is just a number we 
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can move it through anything without any worries): 



^ (3-121) 



D<^°^{k)W^^ = 

l(0) 



(n) " • (3.122) 



Now we note that the operator Dltp is invertible (if it wasn't invertible we wouldn't have been able to derive it in the 
first place - see lecture 7). So if we apply the inverse of D^p (whatever that is - see lecture 7) to the above and we get, 

n^'^fcA = 0. (3.123) 

This is telling us that 11''''' must be proportional to the transverse projection operator P^^ because k\ kills it: k\Pj^ = 
k\ ^5^^ — = k^ — k^ = 0. It is this property which is making the radiative corrections to the photon prop- 

agator cancel out {k^D^^} (fc) liP^D^^^ (k) = 0^ . This is why the photon does not acquire a mass through radiative 
corrections! For example in (/)^ theory the self interaction radiative corrections modify the free <j) propagator such that 
it acquires a mass: 



+, 00 , ^ ^ 



2 .2 

p -A 



Now for the next Ward identity. Consider now H = r]{x)'ip {y) ip (z), 



SeH = {^^'Ap (x)) -0 (y) tp {z) + i] (x) {-ieip {y) euj {y)) -0 (z) + 77 {x) %l> (y) ieeuj (z) -0 {z) 
= e I i {d'^Ap, (x)) ijj (y) tjj (z) - ier] (x) ip (y) oj (z) tp (z) + ier] (x) ip (y) oj (y) ip (z) | 



(3.124) 



Applying our Ward identity {5^11) = we get, 

ei(r (a^^^ (x)) ^ (y) V (z)) - e{Tier] (x) ^ (y) uj (z) (z)) + e{Tter^ (x) i> (y) u (y) ^ {z)) = 0. (3.125) 

We can divide the e away, it's just a (Grassmann) number. Recall that in QED the ghosts are not coupled to any of the 
physical fields in any gauge, this means that the last two vacuum expectation values in the above simplify, 

(Tiei] (x) tp (y) u! (z) tp {z)) = ie{Tr] (x) ui {z)){T-{p (y) tp (z)) (3 126) 
{Tier] (x) ip {y) u {y) (z)) = ie{Tr]{x)u{y)){Tij{y)ip{z)) 

and our Ward identity becomes, 

^ ^<T {d^A^ (x)) ^ (y) V (z)) = ie{Tr] (x) lo {z)){Ti, (y) ^ {z)) - ie{Tr] (x) lj {y)){Ti^ {y) ^P {z)) (3.127) 
This Ward identity will give us a relationship between the fermion 2-point function and the 3-point function: 
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& 




As before let's start by Fourier transforming these vacuum expectation values. We denote the full fermion 2-point 
function (the full fermion propagator) Sp {p) and the 3-point function {k;p,p') (the vertex function). 



(3.128) 



{A^{x){l^{y)^l^{z))= [ 5^4 e^'-^+P y+P'-^) x {27r)U^ {k + p + p') x T^^ {k;p,p') (3.129) 

J (ztt) (ztt) (ztt) 

If we now go ahead and insert these into the Ward identity 3.127 along with the expression for the ghost propagators, 
iv^) = we find, 

/ 7075^75^e<'=-+^-''+^'-) {27rf 5^ {k+p + p') ^F'' {k;p,p') 



{2ttY (27r) 

+ if, r d^k ik.jx-v) i r jlPf,ip.{y-z)c ( \ 



- ie f ^ 



(2j)* 

d_k^ik.{x—z) i 



^J^e^p.iy-z)S,{p) 



= 



(3.130) 



Next we rearrange 



gik.{x-y)gip.{y-z) _ gik.(x-z) ^ik.{z-y) ^ip.(y-z) 
_ ^ik.{x-z)^i{p-k).(y-z) 



and make a shift in the variables of integration to p' = p — fc to yield 

i 

> ' 



'^T^'{k;p,-k-p)+ie—!—SF{p + k)-ie ' 



fc^ + ie 



Sf ip) = 0. 



(3.131) 



(3.132) 




.D^^SF{p)txSF{p + k). 



(3.133) 



Therefore the Ward identity is giving us the result that. 



'k-'+ie^ 



From before we found that ^k^D^"-^ = 



(3.134) 



k^fx - ieSp^ {p) + ieSp^ {p+k)= 0. 



(3.135) 
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Let's see if this makes sense at by plugging in the appropriate tree level values, we have tx = iejx, Sp^ = ^^j^ and 
so we must have, 

ie ]k — re — m) + -e ( — m) = (3.136) 

according to the Ward identity. It is trivial to add up the terms on the left hand side and see that this is in fact the 

case. Beyond tree level we have, Y\ = ie-fx + T'^ where F' is representing the loop corrections to the tree level vertex 
function and for the two point function, Sp^ (p) = ^^j^ + S (p) where S (p) represents the loop corrections to the 
fermion propagator (the self energy). Plugging these into the Ward identity above gives the following relation: 



k^t'x + eY,{p + k)- eS {p) = 0. 



(3.137) 



Chapter 4 



Renormalization and QED. 



4.1 One Loop Correction to the Photon Propagator. 

We now discuss loop corrections and renormalization in QED. Firstly we discuss the photon two-point function: 



which is superficially divergent' (4 — Internal Loop Fermions — 2 x Internal Loop Bosons) = 2 > 0. From 
last lecture we had that the Ward identity gave us that the self energy 11'"' cx P^'^ll (fc^). It is worth noting that 
the BRS symmetry which we used to derive this (and other Ward identities) is independent of the number of space- 
time dimensions. This implies that one can use dimensional regularization techniques to treat ultraviolet divergences. 
Before proceeding with the calculation of the full two point function in dimensional regularization let us first set up 
some conventions. Our gauge fixed QED action in n dimensions is, 

S = j dTx - ^F^uF^"" + + e^j M - mxp^P + ^ + ry^^a^w (4.1) 

We must have that the action is dimensionless, this means the quantities inside the action must have the following 
dimensions (we use [X] to denote the dimensions of X): 

n 

n-2 
2 

^ (4.2) 

= ^-2 
n-(f-2) = 2-f 

We will ultimately want to look at the limit n ^ 4 so (as usual in dimensional regularization) we write, n = 4 — 2e, 
where e is some small positive quantity. Note that e is dimensionless, a dimensionless coupUng constant is a prerequisite 



'The superficial degree is an exercise in counting the number of powers of the loop momentum (fc) appearing in the amplitude for a given process. 
In QED we have a d^fe integration measure and a ^ for every internal fermion and -p- for every internal boson. So in QED we have that a diagram 
has superficial degree of divergence £) if D = 4 — Internal Loop Fermions — Internal Loop Bosons > 



[rf"a;] 
[A] 

m 

[e] 
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of a renormalizable theory. In defining a renormalized theory it is necessary to introduce a mass scale n, known as 
the t'Hooft scaling parameter, which will let us keep e dimensionless through dimensional regularization, we write 
e' = e/i^~^ = e/i^. The /i essentially compensates for the changing of the dimensions of quantities we will calculate 
due to the changing of dimensions of the integration measure, d^x. Introducing the t'Hooft scale parameter changes 
our Feynman rules slightly. 

• Propagators No change, we had calculated free propagators which never depended on e in the first place. 

• "^^^^^^^^"X^ Vertices, iej^^ — » iefi'^j'^. 

The following rules also depended on the dimensionality of space-time and change as follows. 
. /d^x^/d"x,/ (1^^/(0.. 

• g^^Qi^n = n 

• Tr [/] = 4 ^ Tr [I] = n 

• 7''7"7^ = -27" ^ 7'*7"7m = -2 (1 - e) Y 

• -^f^Yl^lti = ^g^" l^YYlti = ^9^" - 2e7^7<^ 

• YYYYliJ, = — 27'^7'^7^ — > YYYY^fj. = —2'y'^YY + 2e'y''YY 

• Traces of 7 matrices — > Not changed unless the trace contains a 7^ which depends on the space-time dimension- 
aUty. 

Now let's get on with our analysis of the photon two point function. Using our (modified) Feyrmian rules we can write 
down the amplitude. 



k+p 



iW {k) = - 




(4.3) 



Counting the powers of the loop momentum p, in the numerator we get n + 2 (n from the integration measure and 
two inside the trace), the denominator clearly has four so the amplitude looks badly divergent in n = 4 but appears to 
exist/be well defined for n < 2. It turns out however that this divergence is an illusion, the Ward identity, 

W"" (fc) = {k^g^^ - k^k,) n (4.4) 

kills the divergence above. We just discussed above that IT'"' has dimensions of momentum n — 2, soinn = 4we have 
[n''''] = 2 in which case the above equation tells us that [11] = 0. 



4.1. One Loop Correction to the Photon Propagator. 



50 



First we use the fact that the trace of an odd number of gamma matrices vanishes to halve the number of terms. 
Finally we use the identities Tr [^i^-y'^-fP^] = 4 {g'^^gP" - g'^^g"" + g'^'^g"'') and Tr [YY] = ^9'^" to give us a 
sum products of n- vectors. 

= Tr [7^ it + l^" iri"" i> + m'^-f'^Y] 

= 4 (kt^p" - gf'^k.p + fc^p^) + 4 (p'^p^ - ^^"'p^ _^ pt^pu^ ^ 4m'^gi^'' 
= 'ik'^p" + Ak^pt" - Ag'^'k.p + Sp'^p" - 4 (p2 _ m^) g^"" 

(4.5) 

The next thing to do is rewrite the denominator in terms of integrals over Feynman parameters, this will eventually let 
us integrate out our loop momentum. To do this we need the identity 

— — ^— = / dxidx2...dxn S {Exi - 1) -—-—^^^—^^ — ^. (4.6) 

AiA2...An Jq [X1A1+X2A2 + ... + XnAn\ 

In our case we have Ai = {{k + p)^ — w?^ and A2 = {p 



dxdy 5{x + y-l) 



1 



1 

7 



(4.7) 



(x((fe+p)^-m2) + (l-x)(p2_m2)) 

Now we try to complete the square in p in the denominator, 



a: ^(/c - TO^^ + (1 - x) (p^ - m^) = xfc^ + xp^ + 2a;A;.p- xm^ - - xp^ + : 

= + 2xk.p + xfc^ - 

= (p + xkf - a;2fc2 + xk^ _ ^2 (4.8) 

= (p + xk)^ + a;fc2 (1 — x) — vn? 

= p2 _ ^2:^2 _ _|_ 

In the last line we have defined p = p + xk, this amounts to a simple constant shift in the integration variable. Seeing 
as the shift is constant we have that the integration measure is invariant / d"p = / d"p and so we can equivalently 
integrate over p, which makes life easier. The photon self energy can therefore be written, 

11 {K) - le H j (27r)" ((fc+p)2-m2-)(p2_„2) 

_ ■ 2 2eflj r cTp 4fe''p''-4xfe''fe''+4/s''p''-4a;/s''fe''-4g^"'fe.p+4a:g'"'fc^ + 8(j5-afe)''(p-x/s)''-4((p-xfc)^-m^)g'"' 

- »e /X Jo da; J {p'^-{xk'^{x-i)+m'^)f 

(4.9) 

As the integral is symmetric in p we can drop terms in the numerator which are odd powers of p. 



n'^-(fc) = ^eV'Jodxj(2^ (^2_(,,,.2(,,_t,+„,2„2 ^ '— 

_ 4- 2 2£ flj r cTp -g'"'p^-2x(l-a:)fc''fc''+2py+g''''(xfc^(l-x)W) 

- M Jo ''^J piy^ (p2_(xfe2(3;_l)+„2))2 

The synnmetry of the integral also then allows one to replace p^p" — > ^p^g^" in the numerator (terms off diagonal in 
H and f look linear in their components with respect to whatever component of p is being integrated over). Defining 
A = xk^ (a: — 1) + we have, 

n-(.)=4^e2,- f\x (4.11) 
Jo J (27r)" (p2_A)2 (p2-A)' 
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These n-dimensional integrals are now in a form that is readily found in tables of integrals,^ 

(4.12) 



r d"i 1 _ (-!)"'» r(m-f ) , I y 

J (27r)" (JS-A)"* ~ (4^)"/2 T(m) VA/ 



J (27r)" (P-A)"» ~ (4^)"/^ 2 r(m) (a^ 

Applying these identities to the two terms in our energy we get (remember in our convention n = 4 — 2e, e = 2 — |^), 

p rf^p 2x(l-x)kt'k" -g''" {m^+xk'^{l-x)) _ i{2x{l-x)k^'k" -gi^" [m'^ +xk^ {l-x))) Y-(e) / i Y 
J (2^Fr (p2_A)^ ~ (47r)^-= r(2) IaJ 

r d"p (|:-i)g''V _ (^2 -i « / 1 Ni-f r4 ni 

J (27r)" (v'^-A\^ ~ \n ^} ^ MttI^/^ 2 r(2) Ia^ V*-'^^) 



(27r)" (p2-A)^ ~ Vn ^ (47r)"/2 2 r(2) 

(47r)"/2 r(2) 



A special property of the Euler gamma function is that zT (z) = T {z + 1) i.e. (l — f ) T (l — -l) = F (2 — |-) . If we 
use this and also the fact that F (2) = 1 second integral above becomes, 

= #.r(e)(i)^A 

Inserting these integrals into 11'^" we have, 

W (fc) = -4e V (47r)'"^ ^ dx {~g'"'A -2x{l- x) ki^k" + g^" (m^ + xk^ (1 - x))) F (e) . (4.15) 

Substituting in for A, 

(fc) = dx (-5'"' (a;fc2 (x - 1) + m^) - 2a; (1 - x) /c'^fc" + gf"" (m^ + xA:^ (1 - x))) F (e) (^)' 

= ^ dx 2x (1 - x) - fc^r) F (e) (^) ' 

= -e (.9- - ^) (^^ /o dx 2x (1 - X) F (e) (^) ^) 

(4.16) 

The photon's self energy is proportional to the transverse projection operator just like the Ward identity from the 
last lecture said it was! Had we believed the magical Ward identity from lecture 8 we could have accelerated the 
calculation dropping all terms not proportional to P}^" = g^^ — The above E^'^ (k) contains divergences as 

n ^ 4 i.e. as e — > 0, we shall now take a closer look at these divergences. For small e, F (e) is approximately given by, 
F (e) « 7 — 7 + O (e) where 7 is Euler's number (w 0.577). We will also be using the following A"^ = exp (e log A) w 
1 + e log A and hence we also approximate F (e) by F (e) « ^ + log A — 7. Taking the Umit n — > 4 i.e. using 
these n we have, 



n(fc2) = lim,^o^/odxx(l-x)F(e)(^)' 

= lim,^o ^ Jo dx X (1 - x) (i + log (^) + log (An) - 7) 

= lim-o f Jo dx x (1 - X) (i + log + log {An) - 7) (4.17) 

= linie^o f Jo^ dx X (1 - X) log + f (i + log (47r) - 7) [^x' - h'] J 

= lim^^o f (Hi + log (4^) - 7) + /o da^ ^ (1 - or) log 



^Note these integrals are in Minkowski space. 
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where we have introduced the fine structure constant a which, in natural units, is equal to 

Clearly the above is logarithmically divergent in the limit n ^ 4 dimensions as is indicated by the ^ outside the 
first term. This is bad at a superficial level clearly as it implies that amplitude for our self energy Feynman diagram 
is infinite which would predict infinity for various measurable quantities e.g. the cross section for e+e~ jJ-'^ ■ 
To obtain a sensible result we have to renormaUze some parameters in the theory i.e. we will the Ultra-Violet infinity 
{i.e. those infinities that occur due to the integrand of a momentum integral becoming infinite as the loop momentum 
becomes infinite) by renormaUzing the parameters of the theory. 
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4.2 Resummation of loops. 

Imagine the QED Lagrangian was rewritten, 

S = / d"a; - lFi,„Ff-' + i^j + eZ^^ /Ail) - mi/^V + ^ (5^^'*)^ + 'qd^'d^u 

= /d"x - \Z^Fn,^,,F'^' + iZ^^R jB^jr + e^jZiV^^ Mr (4.18) 
-mRZoipRipR + 2c|z7 

where ^3,^2,^1,^0, Z^,Zg are real constants and the subscript R denotes a renormaUzed field or parameter. This is 
not the best parametrization perhaps, it would be better to have rewritten the action rescaling explicitly the fields and 
the parameters by constants but we shall see that the above is essentially the same. Consider the part relating to the 
fermions coupled to the field. 

Z2 {i-tpR Pt/JR + eR^-ipR Ari/jr - mR^-iijRil)]^ (4.19) 

We could interpret Z2 as a renormalization of the fermion field tp = \/Z2tpR, which would mean that ^ represents 



Z2 

a renormalization of the mass parameter, m = ^m,R. The renormalization of the charge e and the Af^ field are then 
somehow tied up in ^ as indicated by the second term above. We can untangle the renormalization of e and by 
looking at the first term in the action, 

-IZ^Fr^^^F^-" = -\Z^ (a^A^,. - d^AR,^) - aM^) 

= -\Z^ {{d^AR,,) (aMJi) - {d^AR,) - [d^AR.,) {^^^A'i,) + [d^AR,^) (aM^)) 

= -\ {{dM {d'^A-) - {d^A,) {d-A^) - {d,A^) {d^A^) + {d,A^) {d-A^)) 

— _ip PA"' 

— i^lJ-v^ 

(4.20) 

where we have interpreted the re-parametrization as a renormalization of the A^^ field as A^^ ~ ^/Z^A'^. Consequently, 
from the term coupling the photon and the fermions eZi'il) /Atp we can deduce that the electric charge is being 
renormalized too. In terms of the renormalized quantities this term will equal CRtpR ^RtpR. Hence, 



eip /A'll) = crZiiPr ^RipR 



Z2 V 

^ Z2\/ 



eR^^i)R /Ar-^Ijr . (4.21) 



Finally for the gauge fixing term we have 



_ 1 ^AM /I \2 



id^^A^y . (4.22) 



What have we done? We have multiplied all the terms in the Lagrangian by six real numbers and reinterpreted this 
in terms of a renormalization of the fermion fields, the photon field, the electric charge and the mass. These Z\ are 
where we hide the infinities that come up in the loop diagrams. We interpret S above as being the bare action with hare 
fields and parameters (those quantities with no subscript K). This raises some questions. Does multiplying all these 
constants into the Lagrangian ad hoc break the original gauge invariance (neglecting the gauge fixing term of course)? 
At face value gauge invariance is completely broken, however we have not specified the values of these numbers. We 
shall later essentially demand gauge invariance of this Lagrangian by imposing the Ward identities which will result in 
relations between the Z' s such that the gauge invariance of the Lagrangian is restored. What about the gauge fixing 
term? We do not need to worry about the gauge fixing term, though it is influenced by quantum corrections, as all of our 
S-Matrix elements are gauge independent anyway i.e. we can effectively set Z(^ = \. What about the ghosts? Again we 



4.2. Resummation of loops. 



54 



don't have to worry about what happens to the ghosts as trivially in QED the ghosts are decoupled from the rest of the 
theory and go around as non-interacting (unphysical) complex scalars i.e. there is no renormalization of ghosts as they 
have no interactions Zq = 1. Finally we might worry that we broke BRS invariance? The answer to this is analogous 
to the answer about gauge invariance, i.e. we will impose the Ward identities on the one loop divergent diagrams and 
this is tantamount to demanding BRS invariance as that is where the Ward identities came from. Hopefully the rest of 
this section will give weight to these answers. 

Let's try and work out what our one loop calculation would give using the renormalized Lagrangian. The renor- 
malization should affect the form of the Feynman rules. We derived the photon propagator at the start of lecture 7. To 
do that we rewrote the terms relating to the free photon in the form A^QfjuA", then the propagator was found to be the 
inverse of the differential operator Q^^, in Fourier space. We want to know how our renormalization coefficients affect 
the form of the propagator. Previously our starting point was therefore, 

jd^x-^ (d^A (x)^ - d.A {x)^) {d^A-^ (x) - d-^A^ i^)) ' ^ (^^^'^ i^)f (4-23) 
which has now become, 

Jd^x-^ [d^An {x), - d^An (x)^) {d>^A-^, (x) - ^-^A^ (x)) - (9^^^ {x)f . (4.24) 

We perform the same steps as before to get the propagator. First we put the free photon part of the Lagrangian in the 
form At'Q^^A'' 



I 



d'^x - ^ (^2d^AH,, (x) d'^A'i, {x) - 2d,AR,^ {x) d'^A'i, {x) + ^^d^^R i^) d^A^R (^)) , (4-25) 



and integrate by parts 

(-An ,, ix) a.a^A-^o ix) + Ab_„. ix) d'^d^'An ,. ix) - ^A>i (x) d.A.A'i, {x)] 

(4.26) 



= / d^a; - f {-Ar^, {x) 9^9^ (x) + Ar,^ (x) 9^5'^ A^,, (x) - ^A"^ (x) 9^9, A- (x) 
= Jd'x-f [Ar^, (x) [-g^^di^d, + (l - ^) d^d") Ar,, {: 



Thus once again we find we need to invert the differential operator between the two photon we fields. We proceed 
exactly as before, we need to solve A^^ (fc) Z3 {g^'^'k'^ ~ i}~ CKzl) ^'^^'') = for ^^^^ (^)- 



Where we have defined the projection operators, P^'' = g^'' — ^-j^ and P'^^ = ^-j^, just as before. Consequently, 
as before we have PtPt = Pt, PlPl = Pl and PlPt = = PtPl so Pt and Pl again project out orthogonal 
subspaces. We assume that can be written as a linear combination of these two operators: 

g-i PM = _^ ^pPt^ (4 28) 

Naturally one expects QQ~^ will give the identity or equivalently Pt + Pl'- 

Q ^"Qp.Q-' ^'^ = -Z,e [p^" + ^P^") gp. {AP^ + BP^ 



(4.29) 



-z,e(^AP^^ + ^p^^) = p^^ + p, 



^=-Fk ^ = -|5f (4.30) 
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Free Propagator^^" = D^^ (fc^) 

= -I^APt^' + CrZcPl'') (4.31) 

If we define the full propagator as the sum of all these one-loop corrections, i.e. the free propagator plus the free 
propagator with a loop plus the free propagator with two separate loops plus... ad infinitum we have the Dyson re- 
summed propagator: 

Dap = DoAp (fc')+DoAM {k^) [k') Do.p {k^)+Doxp (fc") iW^-" {P) Do., (fc') ^n-^" (fc^) D^c^p {k')+... 

(4.32) 

This can be expressed in a much more compact form by noting it is of the form 

A{1- By^ C = A{l + B + B^ + B^...) C, 

this is known as resunmiation. To do this we first simphfy the terms of the form Dq" (fc^) iUva (k^) ■ Substituting in 
we have, 

D^'^ in,„ (fc2) = ^ {P-- + (nZcPEnPTuak^ll{k^) 
So the sum for the fuU photon propagator is. 



(4.33) 



- P^J Do„p + (^P^A^Oop (fc') + ^$^P-?A^Oap {k^) + 



{6", - P^,) Do„p [P) + (^1 + ^ + ' + ... j P^A^^Oap (fc^) 

(<5? - ^^^?a) Docp (k^) + -^t^^Qqp (fc^) 

('5? - ^^^a) (-Fk (P^-f + CRZcPLap)) + (-p^ {PTap + CnZiPLap)) 

(^T-ap + (.RZcPhap) + (-Prap + Cii-^C-Piap)) - ^cSh^F)) 

/s2(Z3+n(fe2)) fe223 

-» , £«£i(£3+H(fc!)) p 



(4.34) 



where, 

+ n [e) . Z3 + Inn ^ (i (J + log (4.) ^ ,) + . (1 - .) log ^ ^.^^ _ ,) )) ■ (4-35) 

We can choose so as to cancel the divergences in 11 (fc^) in many different ways, these basically constitute different 
renormalization schemes. We firstly have the condition that this be finite, but how should we fix the finite parts of 
Z^l We could insist that the residue of the propagator was one at its pole k^ = 0. Note that fortunately the radiative 
correction has not shifted the position of the pole - the photon is still massless. Such a procedure is called mass 
shell renormalization. It is a special case of momentum subtraction renormalization. In the momentum subtraction 
renormalization scheme we demand that there be no corrections to the propagator at fc^ = (choosing = gives 
one). Another approach is to fix just so it cancels i, this is called minimal subtraction (MS), or we can choose 
such that it cancels the 7 — 7 + log 47r which is known as the modified momentum subtraction or MS scheme. In these 
renormaUzation schemes the scale jx is left over, it's numerical value may be chosen at will - in relationships between 
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physical measurements / processes it drops out. In MS, 

Zs = -^(i-7 + log47r) 



(4.36) 



We should also choose = Zs to render the longitudinal part of the fuU propagator finite (see last line of 4.34) ! Then, 
letting 

n(t') = f|d».(l-.)lo,( ^,^^^:,^_^, ) ,4.37, 
we have for the one loop corrected photon propagator: 

Note that the pole has not shifted so the photon has not acquired a mass through these radiative corrections but instead 
the gauge has shifted! The fact that it is attached to the longitudinal part of the propagator is telhng us that we should 
not expect contributions to matrix elements from the longitudinal part of the resummed propagator. Also for fc^ > 4m^, 
n(fc^) has an imaginary piece, this is associated with the fact that for > 4m^ there is enough energy to create real 
fermion-antifermion pairs in the loop. 
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4.3 The Electron Self Energy and the Vertex Function. 

In this section we will return to the Lagrangian and Feynman rules of the bare Lagrangian i.e. everything is bare, there 
are no Z's. At the end of the section we will shift back to the renormaUzed parametrization of 4.2, look out for that. 
We could proceed to calculate the one loop corrections to the fermion self energy: 



p p-k p 



and the vertex function. 




In terms of the self energy the fuU fermion propagator satisfies: 



If we use the Ward identity 



we find that 



Sy = Sp^ + iHip) (4.39) 
Sp^ = -i (i. - m + S (p)) 

A;^f> {k,p) = -eSp {p+k)-^ + eSp (p)"' (4.40) 

fc^I^(fc,p) = ze(S(p + fc)-E(p)) (4.41) 



and remember that this is true in any dimension i.e. for arbitrary e. Bearing this in mind let us look a httle more closely 
at the structure of the divergent diagrams. 
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/ 



A 7"^^ 7^ . , , ,\ 7'% (4.42) 



where P^x is the photon tensor structure (for whichever gauge) and is dimensionless. Rewriting the integral, 



we see that, only when both /q factors are taken in the numerator is it divergent. Doing the integral we know that 
/q... /q = qa-.-qrl" ■■■Y ■ Note that 7'*'7''7'^ = {2-n)Y /q... /q = ^-^^ ...Y ■ So the structure of the 
divergent term must be simply 7** ! Other pieces proportional to p^^,k^, rwy^ etc must be finite ! 

= ieY (j- + h (p^) j + other finite terms. (4.44) 

Now let us turn our attention to the fermion self energy: 

Now this integral is apparently linearly divergent; however we know that this must actually integrate to zero leaving us 
with logarithmic divergences. We will find that, 

S ^^i* (7 + 5 (/)) + m (I + / (4.46) 
where g and / are regular functions. Returning to the Ward identity and substituting in we find that 

iekf,{Y{^+h{p^))+ finite terms) = +ie{{^+ /:) (7 + .9 ((p + fc)')) + m (f + / ((p + fc)')) 

- ^>{i+9{p'))- m{^+f{p^))} 
=> ie /c^ + finite terms = +ie /c^j + finite terms 

(4.47) 

Since this is true for arbitrary e we conclude that C = +A. That is to say, the divergent term in the proper vertex is 
related to the divergent term in the wave function renormalization of the fermion self energy. 

1 

Now we flip back to the renormalized parametrization 4.18 multiply the (bare) fermion fields by (in notation of 
section 4.2 Z2 = Zj), 

^ = ZJtPr (4.48) 



and write, 



then 



must be finite. In MS 



m = —mR=mR + 6m, (4.49) 
A B \ 

Sp^ = -i[{i)-mR - 5m) Zf + ^- + mR— + ... \ (4.50) 

A 

Zf = l (4.51) 



e 



then ZfmR + 6m = mR + l) i-e. 



-5m = ( 1 - - ) mfl - ( 1 + — ) Wfl = -mR ^^^ . (4.52) 



Now for the interaction term we have. 



e^R MnZf^fzl (4.53) 
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finally we multiply e by (= „ •^^^^ in notation of section 4.2) to give, 



eRi>R MrZ^Z^^Z^ (4.54) 

so that ^67'* [ZjZ^^'Z^ + ? + finite) is finite. However, C = A = e — eZf so to cancel the infinities here we 
require, 

, / ^ (4.55) 

= ^ = 

Z^Ze = 1 (4.56) 
In the notation of the last section ( Zg = 1 I this result is, 

Zi = Z2. (4.57) 
To summarize we see that the Ward identities mean that: 

• the terms multiplying ^/i /)ip and e'tp /Aip in £ are the same 4.57. 

• the charge renormalization {Z^) depends only on the photon wave function renormaUzation [y/Z^j . 

From the discussions of sections 4.2 and 4.3 we have trivially Zq = 1 and not trivially Zi = Z2, Zt^ = Z^. Substituting 
these expressions into the bare action we have, 

S = y d"a; - -^Z:,Fr^^,F^^ + Z.^Jr {i p + e ji A) /Pr - niRZoi^R^R + ^ (9^^^)' + vd^d^u; (4.58) 

If we now write the action in terms of a renormahzed action plus a counter-term Lagrangian viz S = Sr + Sct, 
we have. 



Sr = J d"a; - iF,i,;..Fr + fiJ (M + e A) /pR - mR^^Ri^R + ^ (9^ A^)' + nd^d^uj 
Sct = Jd^x - ISsFr^^^F^"" + S^i^R {i jd + e /A) /pR - SmiiR^R 

5i = Zi — 1 Sm = Zq — I 
^3 = Z3-I 



(4.59) 



(4.60) 



Note the ghost plus gauge fixing sector here is the same as in chapter 3. In that chapter we saw that the basic part 
of the Lagrangian (kinetic and mass terms for physical particles and their interaction) is BRS invariant because gauge 
transformations are contained within BRS transformations (albeit with a fancy function comprised of Grassmann num- 
ber times Grassmann function) and that the gauge fixing term plus ghost sector are together also BRS invariant. The 
gauge fixing and ghost sector is therefore BRS invariant, as before. The photon kinetic term was gauge invariant 
(and so BRS invariant) and so is the renormahzed version in Sr but now under (renormalized) transformations of the 
same form as the original ones but with unrenormalized quantities replaced by renormalized ones. The same goes for 
V'fl jS + ^R A) i'R mRtjjRtpR, they continue to be separately gauge invariant under the gauge transformations 
of the same form as before but with unrenormahzed quantities replaced by renormalized ones. The renormahzed action 
is invariant under the original gauge transformations with unrenormalized quantities replaced by renormalized ones, 
the counter-term Lagrangian is invariant under renormalized gauge transformations. 

We now recap what has gone on and add some words of warning. In the case of a theory with a global symmetry 
it is easy to show that the counter-term Lagrangian is also invariant under the symmetry (see e.g. [5, 8, 13]). One might 
therefore expect that the counter-term Lagrangian of gauge theories would behave hkewise. This isn't true though, 
gauge fixing broke gauge synnmetry long before any talk of renormahzation and showed us that the true synnmetry of 
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gauge theory is in fact the BRS invariance related to the gauge group of the invariance of the un-fixed Lagrangian. 
So is the counter- term Lagrangian supposed to be BRS invariant? In the case of a general gauge theory i.e. non- 
Abelian this is not the case either, the obvious symmetries that we can see in the counter-term and renormalized actions 
are unique due to the Abelian nature of QED. The coefficients of gauge variant counter-terms vanish in QED. In the 
general case the gauge variant counter-terms can arise, they are such that the renormalized Lagrangian is invariant 
under renormalized BRS transformations. For more on the general case see [3,4]. We would also like to emphasize 
that what has gone before is not a proof of the renormahzability of QED. Here by renormalization we mean not just 
in the power counting sense of dimensionless couplings but that renormalization and gauge invariance are compatible, 
that the renormalized action is gauge invariant under some representation of the gauge group which the bare action is 
invariant under. Our treatment has been merely exploratory, a tour of renormahzability in the sense just described. We 
arrive at a gauge invariant renormalized Lagrangian by inferring Zi = Zi which was in turn derived from demanding 
that the action including the renormalization constants Zi obey the Ward identities, which was equivalent to demanding 
the renormalized action be gauge invariant. A nicer (longer) thing to do would have been to calculate explicitly, with a 
gauge invariant regularization, the graphs above and show the infinities cancel naturally rather than demanding it occur. 
Formally the proof of renormahzability is an inductive one. It is possible to expand the (effective) bare action in powers 
of h which is equivalent to expanding the action in loops, the power of h in the expansion corresponds the number of 
loops associated with that order. The induction proof requires that one show that the loop expansion at n loops and at 
n-\-\ loops obeys the Slavnov-Taylor identities^. In other words the proof of "renormahzability" amounts to showing 
that the (BRS) symmetry, expressed by the Slavnov-Taylor identities, exists at each order of the loop expansion. In fact 
the BRS transformations are not necessary, it is possible to obtain Ward / Slavnov-Taylor identities without it but the 
BRS machinery greatly simplifies them which in turn makes the problem of proof of renormahzability of gauge theory 
a tractable one. We recommend the reader to explore this technology further in the literature of some of its founding 
fathers J.C.Taylor, B.W.Lee, J.Zinn- Justin and G.'tHooft [2, 3, 6, 13]. 



^Slavnov-Taylor identities are what Ward identities are called in the context of non-Abelian gauge theories, they are derived by analogy to 
what we have been doing so far. Generally "Ward-Slavnov-Taylor identities", "Slavnov-Taylor identities" and "Ward identities" are somewhat 
interchangeable names. 



Chapter 5 

Anomalies. 



5.1 Chiral Symmetry^ . 



We just observed that there is no radiative correction to fermion masses if the mass parameter m ^ 0. Thus a massless 
fermion will remain massless to all orders in perturbation theory, behind this lies chiral symmetry. Chiral symmetry is 
an internal symmetry. We shall see that a massless free theory (in even dimensions) possesses U {!) chiral symmetry. 
To begin with we will study chiral symmetry in a slightly modified version of QED, axial electrodynamics: 



S = Jd'^xC 

= fd^x- iF^^Ft^" - IG^.^G'"' +ip{ifi + qy + g /A-f' - m) V + Ghosts + Gauge Fixing 



where F^i, = d^Vv — d^V^ and G^^ = d^A^ — di^A/^. First consider local vector gauge transformations Uy (1) 

=> Vi ^ ^g-i9a(a;) (52) 

^ V^ + d^a{x) 

These are essentially the regular U {!) gauge transformations and they only act on what is basically the QED La- 
grangian i.e. we know that this will leave things invariant by inspection barring the term ■tpi'y^igAn'y^ip which is also 
trivially invariant: 

We also know that by using the classical equations of motion the symmetry transformations of the fermions (i.e. only 
the first two transformations above) result in the conserved vector current jy = tpj'^ip so-called because it transforms 
as a 4-vector under Lorentz transformations. So we have local Uv (1) gauge invariance, does the Lagrangian also have 
C/^ (1) gauge invariance i.e. is the Lagrangian invariant under, 

=^ V' ^ ^g+Pis/3(a;)7' (5.4) 

A^, A^ + Qd^P{x) 

? We have generalized the transformations with real numbers P and Q for reasons we will come to. Let us work 
through the terms in the Lagrangian. Clearly the kinetic term —jF^j^^F^^ of the vector field is unaffected so we 
need not consider it. The term — jG^^G^^ has basically the same composition and transformation as —jF^^F^^^ 



'Unless otherwise stated the terms 'chiral symmetry' and 'axial symmetry' are interchangeable in these notes. 
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which is itself invariant so we can also add it to our list of invariants and we need only concern ourselves with the 
•ijj{...)tp part. To do this we wiU need the following trivial identities, 



^^e^n-h' = 7m (l + if (^) 7' + I {^f {x) 7^)' + = (l " (^) 7^ + ^ {if i^) 7')' + 7^ = e'^^^^^^^^ 



(5.5) 

Q^e'fi-h' = ^Q^j ^^^^ _ (5 7) 

Now we use these identities in considering the rest of the transformation, 

^{iP + qy + g /A-f' - m) V ^ ^gPis/3(x)7' ^ + ^ ^ + ^ j^g{g i^x)) 7^ - m) e^'sis{x)^^^ 
+ (g /l^ + 5 /Ay' + (x) 7^ - m) e^^^akx)-^'' ^ 

^^QPigf3(xh' f^^-Pig0{xh'^5(^gp ^jS (x)) ^ + le-^'S/Bixh' 

+ e-P'3l}{x)j' {q y + g /Ay' + gQ ^0 {x) 7^) V - me^*»/5(^)T' v) 

= (75 (ffP (a;)) ip + i jd^ + q y^ + g /Aj^ip + {gQ jd/3 (x)) 7^^ - me'^^'sdixh' ^'^ 
= ^l){ijd + qy + g /Ay' - m) 'ip + {\ - e2«9/3(^)7"^ ^ + (^q _ p) g (5^/3 (a;)) ^^/^^s^ 

= Vi (i ^ + g + 5 ^7' - m) V + mVi (1 - e2«s/3(-)T' j V> - 5 (Q - P) /3 (x) + g (Q - P) 9^ (/3 (x) j>^) 

(5.8) 

We have defined the axial current = ^jj-f''j^tp so-called because it transforms as an axial vector under Lorentz 
transformations i.e. when a; — > — x, — > j^. In the action, the last term in the above is a total divergence which 
we can rewrite as a surface term, the surface being at infinity. Assuming that the fields and their first derivatives are 
vanishing at infinity we set J d^a; g{Q — P) (/? (a;) j^) to zero - we did this aU the time at the start of the notes. To 
first order in l3 {x) we have^ 

S^S + j A^xgl3 (x) ((P - Q) d^j'X - 2imPij-r'^ij) . (5.9) 

If we take the standard gauge transformations P = Q = 1 we have that the action is invariant up to J d^x 2im^7^ V> 
so in the Umit of massless fermions the action is invariant under local chiral gauge transformations. Massless fermions 
are a general feature of chiral invariant actions. On the other hand had we just looked at local gauge transformations of 
the fermion fields i.e. P = 1, Q = such that Af, A^ we have from demanding invariance of the action that the 
axial current must obey the following conservation law, 

d^.j'X = 2imi}^^i) (5.10) 

as (3 {x) is an arbitrary function, ip^^tp is a pseudoscalar under parity operations x —x it transforms as tpj^tp 
—'ipj^tp, the four divergence of the axial current is proportional to the mass density of the pseudoscalar. At the classical 
level we can show that this current relation holds by considering what happens with the classical equations of motion 
for ip and ^. HopefuUy getting the equation of motion for tp is trivial, to get the equation of motion for •ip one can 
conjugate that for ip using the fact that 7°, 7^7'' and 7^ are Hermitian (see below) or rewrite the Lagrangian using the 
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product rule and vanishing surface terms so that the derivative in the kinetic term acts instead on ij;} 

[i /9+ y+ /A-f- - m) = 

^ V^7° {-K + V + - rn) (707/') 7O = (5.11) 
^ ^ (-iS^ + F- A75 -m) 7" (7"7f)7° = 
^ (i'^- /l^- /I75 + = 

To derive the classical current equation we use these equations to substitute into the current equation ^ ^ and fill): 



= -iVi ( y+ /I75 - m) 7^^^ + 11/^7^ (- y- /A'f + m)ip 
= y+ /A'X' - m) j^ip + i^j) ( y+ /A'T' + m) y'ii 

= 2im'tp'y^'tjj 



(5.12) 



At the quantum level what we have done is effectively redefine our variables of integration (the fields) by ^ — > 
■^e^l^i'^)!^ and ij: — > e^^^^'^^^ip so that, (omitting spectator fields for brevity), 

= J DtPDtl) exp i^d^xS [i), V] + 5/3 (a;) (S^Ja - '^irmp-^^tl)) + O {(3^) 

(5.13) 

We have assumed (naively) that because. 



Det 



Det 



V 





-i/3(x) 






\ 



ei/3(x) 

e*'^^^) y 



the Jacobian associated with the path integral measure under the transformation is also one hence. 



(5.14) 



(5.15) 



This gives us the {naive) axial vector current Ward identity. It is naive because it is only true at the classical level, in 
our abbreviated path integral description above we neglect the possibiUty that the path integral measure may change 
non-trivially under the transformation of the fields, it may contribute something in the form of a Jacobian. When the 
Jacobian is not one we have an anomaly. Anomalies are of fundamental importance in field theory. 



^dfi indicates that 9^ is acting on all the stuff to the left instead of the right. 
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5.2 The ABJ Anomaly. 

Before plunging into a calculation of the change in the path integral measure we will detour to some phenomenology 
which preceded it. As we said before the anomaly presents itself by making the path integral measure over the f ermions 
transform so as to produce a Jacobian not equal to one. It also presents itself in the calculation of certain S-matrix 
elements, the axial anomaly was first discovered in the calculation of the VVA (vector, vector, axial) and VVP (vector, 
vector, pseudoscalar) triangle diagrams. 




These correspond to the following 3-point functions, 

{Q\Tj^{x)j,{y)jA,x{z)\0) (5.16) 
{0\Tj^{x)j,{y)P{z)\0) (5.17) 

In momentum space these are, 

7;.A (fci, fc2, q) = ij d^xd^yd^z e'krMy-igz ^0 |^_^.^ ^^^j^ ^y^-^^^ (^^1 0^ (518) 

T;.(fci,A;2,g) = i j d*a;dV'^ 6''=^"+'"^™ {Q\Tj^{x)j,{y)P{z)\Q) (5.19) 

with P {z) the pseudoscalar i/; (z) y^^p (z) (probability) density. We are trying to test our naive chiral ward identity. 
The first amplitude contains an axial current and the second contains a pseudoscalar, the Ward identity relates the 
divergence of the axial vector to the pseudoscalar density, so should provide a relation between these two amplitudes. 
To do this we will clearly want to be differentiating the first amplitude so as to get the divergence of the axial vector 
current. This will involve differentiating the time ordered product of operators (x) jv {y) jA,x {z) which is not easy. 

( e{xo- yo) e (yo - ^o) 3^ {x) j" {v) 3\ {z) + e{xo- zo) 9 {zo - yo) F (x) j\ (z) f {y) 

= 0{yo-xo)9{xo-zo),r{y),r{x)j\{z) + 9 [yo - Zo) 6 {zo - xo) j- (y) j\{z) (x) 

[ + 9{zo-xo)eixo-yo)j^{z)j''{x)r{y) + 9 {z^ - yo)9 {yo - x^) j\{z) f {y) j^' {x) 
= [Tj^(x).r(2/) {dxj\{z))] 

- 5{zo-yo)9{xo-yo)f{x)3-{y)j\{z) - 5 {z^ - xo)9 {zo - yo) f {x) jliz) f {y) 

- 6{zo-xo)9{yo-xo)r{y)j^{x)j\{z) - 6 {zo - yo) 9 {zo - xo) {y) j\{z) j'' {x) 
+ S{zo-xo)9{xo-yo)fA{^)j''{x)r{y) + ^ {zo- yo)9 {yo - xo) f^z) f {y) {x) 
+ 5{zo^yo)9{x^-z^)j^^{x)j'i{z)f{y) + 5 {zo - xo)9 {yo - zo) f {y) fA{z) {x) 

= [rj^(x)j-(y) {dxj\{z))] 
+ <5 (zo - yo) {9 (yo - xo) 4 {z) f {y) {x) ~ 9 {xo - yo) {x) f {y) j\ (z)) 
+ 5{zo- xo) {9 {yo - zo) f {y) j\ (^) j'' {x) - 9 {zo - yo) {x) {z) f {y)) 
+ 5{zo- xo) {9 {xo - yo) {z) {x) f {y) -9{yo-xo) f {y) {x) jl {z)) 
+ 5{zo- yo) [9 {xo - zo) j'* {x) j\ {z) f {y) -9{zo- xo) {y) A (^) 3" {x)) 

(5.20) 

In actual fact the thing that we are differentiating is inside an integral (it's inside Tp,yA (ki .k2,q) at the top of the page) 
so when we do the integration the delta functions will set some of the xos, yos and zos equal to each other. We can 
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therefore safely rewrite the 6 functions as if this has already occurred. 

= [ra.^j^^ {x))f{y)j\{z)\ 

+ 5{zo- yo) {e [yo - xo) j\ (^) f (y) {x) -0{xo- yo) (x) ^ iv) j'a (^)) 

+ 5{zo-xo) (e {yo - xo) f (y) (z) (x) -9{xo- yo) (x) (^) f iv)) 

+ S{zo-xo) {e {xo - yo) fA (^) (x) r {y) -0{yo- xo) j" iv) (x) j\ (z)) 

+ 5{zo- yo) [e (xo - yo) F (x) (z) f (y) -6{yo-xo) f (y) 4 (^) J'^ i^)) 
= [T{d,^^j^ ix))riy)j\iz)] 

+ (0 (yo - xo) {6 {zo - yo) [fA (^) , f iv)] ) (■^) + i^o - yo) {x) {6 {zo - yo) [4 (^) , (y)] )) 
+ {9 {yo - Xo) {y) {S {zo - xo) [j1 {z) , 3^ {x)\ ) + ^ (xo - yo) (^o - ^o) \j\ {z) , 3^ {x)\ ) 3" (y)) 

(5.21) 

Now consider the commutators, they are all accompanied by a delta function which gets integrated over in Tp,^ a ( fci , ^2 , g) 
effectively making them equal time commutators. We know the equal time commutation relations (ETCRs) for fermion 
field operators such as those above from canonical quantization: 

{V^„ (.t) , V/3 (2/)} = 

{^iix),4{y)} = . (5.22) 

{ (a;) , V'/3 {y) } = SociiS^ (x - y) 

The a and (3 are Dirac indices. Given these ETCRs consider a commutator of the form, 

{x) Kp^P (^) . ^\ iy) Ks^s {y)] (5.23) 

Where F'^ and A'' can be 7*^ or 7*^7^. Written in terms of Dirac indices the F'' and A" are just numbers (as opposed to 
operators - tp^s and ips) and can be puUed through as such. 

= {^i {x) # {x) {y) i^s {y) - i'l {y) i's {y) € i^) # {x)) K^A^s 

Pull the V'a (x) ipis {x) through to the other side of the first term. 
= {S^ (x - y) Sff^tpl {x) ips {y) - 53 (x - y) 5ocStp\ {y) ipp {x) 

+ i>\{y)i^s{y)^i{x)^0{x)-^\{y)i,s{y)i>i{x)i>p{x))V''^^K';^ (5.24) 
= 53 _ y) (^t (a,) r^^A;^,^^ {y) - i;i {y) A'-^^T'^^i^s {x)) 

+ (V-t (y) {y) v4 (x) V/3 {x) - V4 {y) {y) ^1 {x) i^p {x)) r^^A!;, 

= (53 (x - y) {x) F^^A^.V^ {y) - V'i {y) A^^F^.V^ {x)) 

Given that the above is multiplied by 5 {xq — yo) (x — y) and we integrate over x and y, after one of the integrations 
this becomes, 

V'^a;) [F'',A'']V(a;). (5.25) 

In our particular case we always have a 3^ {x) (= ip^ {x) tp{x)) in the commutator i.e. in terms of the above we 
are always deaUng with 

everything commutes with the identity so all our commutators vanish giving. 



i^) r°/3V'/3 {x) , {y) A-giJs {y) = {x) [r^, A'^] V (a;) with F^ = / the identity. 



a,%.A(fci,fe,g) = »/d4xdV^e*i-+^'=^™ ' ^9 (0|Ti^(a;)>(y)iA,A(^)|0) 



= i J d''a;d''yd^2; e''=i^+''=22'-'9^ 



+ d^{0\Tj^{x)jAy)jA,x{z)\0) J 
+ {0\Tj^{x)j,{y){d^jA,x{z))\0) 



(5.26) 



We can now use the naive Ward identity to rewrite the 4-divergence of the axial current as d^jA,\ {z) = 2imP {z) 
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where P {z) =il) {z) 'y^tjj {z) is the pseudoscalar (probabiUty) density. 

' ^ ^ \+ 2im{0\Tj,ix)jAy)P{z)\0) J 

In fact d^T^i^x (fci ,k2,q) is an integration over a 4-divergence which we can take to vanish using the same trick were 
we write it as an integral over a surface at infinity and assume all the fields and their derivatives vanish at infinity. 
Hence, 

J d^xd^d^z 2m(0 \Tj^ {x)3. {v) P {z)\ 0) - q\0 \Tj^ {x) 3. iy)jA,x {z)\0) (5.28) 
Consequently we get the naive Ward identity for our two three point functions, 

q^T^.X (fci, fca, q) = 2mT^, (fci, fca, q) • (5.29) 

It is also worth pointing out that there are also vector Ward identities associated with the VVA amplitude which are 
derived in exactly the same way as the axial Ward identity. In this case we differentiate with respect to one of the other 
coordinates in the time ordered product, x or y. This gives a time ordered product of the divergence of one of the vector 
currents with the other two currents and a bunch of equal time commutators of the form [j*^, j'*] all of which vanish 
just as they did in calculating the axial Ward identity. Applying the classical conservation of the vector current to the 
time ordered product gives the vector Ward identities for the VVA diagram, 

KTiJ.'^\{ki,k2,q) = 
k^T^^x{ki,k2,q) = ■ 

Hopefully it's clear that for V VV diagrams we would get three of these and for VAA we would get one vector and two 
axial Ward identities. 

Thus we have a Ward identity relating the VVP and VVA amplitudes which was derived assuming that the classical 
chiral Ward identity d^j'^ = 2imP remained true at the quantum level. We can now test this hypothesis by simply 
evaluating the two diagrams and see if the relation above is obeyed. It turns out that the divergence of the axial current 
does not obey the classical relation. The classical relation is modified by an anomaly and the Ward identity above 
relating our three point functions is missing a term A^i, = — 2^^fiiyai3kfk2 on the right, the so-called ABJ anomaly 
(Adler [1969] Bell and Jackiw [1969]). 

Hopefully it is obvious that the Feynman rule for the axial current vertex is just e'y^j^ as opposed to 67'* for the 
usual vector current vertex. With the usual Feyrmian rules we then get the following the ampUtudes, 



T^uX {ki,k2,q) 



Tij,u {ki,k2,q) 



To check the Ward identities we can use the following, 

{A7'}=0={i',7'}-2m7' + 2m7' (5.32) 



/^x-lxTr 



/(l^x-lxTr 



1 ^ i ^ 



-L. ' /yO i ^ L 



(5.31) 



(5.33) 
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gives, 



1 ~1 " 



_ r d^p J, 

+ 2mJjp^x-lxTr 
= 4-;^') + + 2mT^, (fci , , g) 



y(—m, ' 2'^—<~f—m ^ -^^ --- 

i ^5 i i 

2/— m ' ^—</—m 



*7 j/.g^-m^i^y-^i- 

+ ^7 j/-^-to7;*^l:^ 



r7^ 
77^ 



iJ^t"^^^ and -R^y^^ are the so-called rej? ienw*. 



(5.34) 



'\iv'' — if ■^:;^Tr 

,iA2) _ . r J^rr. 



We can 



+ ^^7^7M2/-^2-m7i' 

7 7i^ _m 7/i 

L 7 ^-^-m^lJ. j/_^2_m7i^ 

rewrite i?^'^^^ with {7^, 7^*} and the cyclic property of the trace so as to closer resemble iJ^'i 

1 s 



'-m 7^7i' _m 7m 



(5.35) 



,(A2) 



(27r)' 



+ 



+ 7^ 



1 5 ^ 

^-m^ivl y_^^_m7M 
I /V 1 /V 

T^-mil^l j/-^2-m7f 
51 1 

5 



(5.36) 



So we have the axial current Ward identity if the rest terms, 



{A2) _ 



7^ ^-}f^-m^v 

+ 7^ _m 7ai 'ff^lv 

7 -^2-m 7"^ j^-^i -m 7m 

"I" 7 j/- ^1 - ^2 - m 7m ^2 - m 7:^ 



(5.37) 



cancel each other (they're clearly not both zero). As written above it looks like simply shifting the integration variable 
p by a constant to p + fci in the first term in r\^^^ and to p + A;2 in the second term in r]^"^^ we would have 
R^Mif^ = —-R^y^\ a cancellation of the rest terms and agreement with the axial Ward identity! This is not true however 
as these integrals are divergent, shifting the variable of integration in divergent integrals is not generally trivial and 
there may be a finite difference between the shifted and unshifted values. It is equivalent to saying that the amplitude 
T^vX {ki ,k2,q) depends on whether we use p for our loop momenta or p + constant, naturally one expects this not 
to matter, p is a dummy variable and the hmits of integration are at infinity, however in this way the amplitude is 
ambiguous. We have a similar situation for the vector Ward identities k^T^^x {ki,k2,q) = kl^T^^^x {ki,k2,q) = 0. To 
see this we use the identities 



i»- m 



= 1 - (- i»- m) 



1 



and 



1 



i» - ro i)- /i-m 



/l-m 



(5.38) 
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(2^) 



_ r d p J, 



= -I 



Tr 



-^—ra 
i 



1\T 
1\T 



'^—cf—rn /'^ '{/—J^i—rn 



i 



-g'-m 

«5 



5 



—m 



j/-|«2-m 



Tr 



'T-^'^ '^—(l—'m 
■ 1 5 1 



- ^1\T 



1^—1/1—1/2 — m ''''' J''— ^1 — ™ 



7^; 



+ Tr 



m) ■ 



;;7i' 



(y2) 

Ai/ 



(5.39) 

An analogous calculation gives a similar result for k^T^t.v\ (fci , fc2 , g), in both cases we see that shifting the variable of 
integration in one of the terms by a constant would give a cancellation and the classical vector current Ward identity. 
Like the axial case we cannot trivially do this as the integrals are linearly divergent. We will now calculate the rest 
terms and treat properly the shift of the loop momentum. 

We will first attempt to calculate the rest terms in the axial case starting with the trace. The trace of the first term in 
R^^^^ is 



Tr [7^ (^- Jki + m) 7,. (^- ]kx + m) 7^] 
Tr [75 (i-- JiY- h) 1^ {i>- h) In + "^^7^7.^7^] 
Tr [7^7.7^] = and Tr[{i>- q) 7. {i>- h) 7^7^] 
4iea/3i.^ {p - g)" (p - ki)^ 



= -AieauPiJ. {p - q)" ip - kiY 



(5.40) 



r/iAl) 



(27r)* p^-m^ 



(p-fc2)°P^ 

(p-fe2)^-m2 

(p-fcl)°p'' 

(p—ki)'^—m^ 



d^p 



(p-fcl 



(2^)^(p-fei-fe2)''-m2 



(p-fci)" 

(p-fc,-fc3)°(p-fc^)P 



(5.41) 



{p-k2V- 

It is easy to check that terms of the form eapfj.uP'^P^ are identically equal to zero this is what kills the apparent quadratic 
divergence making the integrals linearly divergent instead. To get a cancellation of the rest terms above we need to 

(Al) (Al) 

shift p — > p — fci in the first term of i?/,,^ ' and p — > p — A;2 in the second term of Rj^v ■ We can write the rest terms as 
functions of the loop momentum, we parametrize the first case as having no shift i.e. 



TfiAl) 

p>{A2) 



^J^{A,Ap)-ki^k2) 

+ J^iA,Ap-ki)-k^^k2) 



(5.42) 



where the As correspond to the first and second terms in i?^"^^"* and R^^"^^ . So if the integrals were convergent such 

,4 ,4 
shifts would have no effects, we would find J j^:^ (p) — fci ^ k2) = J (^^'^ (P ^ ^1) ^ ^ ^2) and 

the rest terms would cancel but the integrals are divergent so this is not the case. To find the anomaly we need to know 

the extent to which Rln, and Rln^ don't cancel i.e. we must determine how the rest terms are affected by shifts in 

the loop momentum. 



A. ^ f ^"P. (p+a-fc2)°(p+a)^ 

'iea/J/iL/J (2^(p2_„2)((p_fe2)2_„2) 



(5.43) 



We can Taylor expand A^^ (p + a) about A^^ (p) to give. 



+ a) - A^^ (p) = oT-d^A^^ (p) + a'^a^a^^A A^^ (p) + 



(5.44) 
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hence, 



= J^Ia dS^'^r (A^, (p) + a^dxA^, (p) + ...) 

where in the last line we have used Gauss' law to rewrite the volume integral / d^p 9^ {o.'^fnv {p)) as a surface inte- 
gral dS'^Grffj.iy (p). The surface ^ is a 4d (hyper)sphere of infinite radius (the volume being integrated over was 
symmetric and infinite). The integration measure dS'^ is an infinitesimal 4-vector defined on A and perpendicular to 
it. Now to evaluate the surface integral (we replace {p — k2)" p^ with — fcfp^ in the numerators of the integrand as 
eafs^..p"p^ = 0). 



■ / dS^ar (A^. ip) + a'dxA,, (p) + ...) = / dS^a^ 

Ja (ztt) Ja 



( \ 

(p2_TO2)((p-fe2)^-m2) 
^ ^ (p2_TO2)((p-fe2)^-m2) 



V + ) 

(5.46) 

As the surface of the sphere we are integrating over is at infinity we can greatly simplify the denominators above by 
taking the limit — > oo. 



-1-4 / dS-^a, (A^. (p) + a^dxA^, (p) + ...) = J^I^iH^ f dS^a^ + a'd, (4) + -) (5-47) 
(ztt) Ja (27r) Ja \P \P J J 

In three dimensions the integration measure is easily found to be rfS*' = |p| sin9 dOdcf), where = p is a vector 
from the centre to some point on the surface of the sphere (i.e. it is radial) and 6, 4> are the usual spherical polar and 
azimuthal angles. Analogous to this in four dimensions the integration measure is proportional to p"^ |p|^ consequently 
the second term above does not contribute as it has fewer and fewer powers of p which makes that integrand tend 
to zero all over the surface of integration (|p| oo). The same is true for the other higher derivative terms in the 
Taylor expansion. However we can see that the first term contains an equal number of powers of p in the numerator 
and denominator in the integrand which making it non-vanishing, we will evaluate this below. Had the integral been 
quadratically divergent the second term would have contributed, cubic divergences would have made the third term 
contribute etc. Following this reasoning we find that only in the case of integrals which are convergent or at most 
logarithmically divergent does a shift in the integration variable not contribute any such terms and the final integration 
is unaffected. We can also generaUze our 3-d spherical surface element to a 4-d one, by analogy it will be of the form 
dS"^ = p^p^f {9, 4>) dOdcjxi'tp where 9, (p and are the spherical polar angles of Minkowski space. Also the 4-vector 
a'^ is a constant and so it can be brought outside the integral. 

=^ j^lAdS-ar{A^.{p) + a^dxA^^{p) + ...) = /^P^ \pf f (0, <t>) a.g^ dgd#V 

Needless to say an honest evaluation of this integral is highly tedious. Instead of showing the 16 integrations we sketch 
how they are done. First to properly determine the integration measure / [9, ^p) d9d(j)dtp one has to evaluate / {9, 0), 
this is the infinitesimal surface area element. It is equivalent to the determinant of the metric (of the sphere). In 3-d the 
metric of a sphere of radius R is, 

f gee 9e,\f \ ^^^^^ 

\ 94>e 94>ci> ) \ Q R^sin^e J 

To get the metric all you have to do is take a point on the sphere parametrized in spherical polars and make infinitesimal 
rotations in all the angles then compute ds^, the infinitesimal distance between the new and old points. The root of 



the determinant of the metric above is V-R^ x R'^sin'^9 which gives the well known R^sin9 d9d(t) as the surface are 
element. PictoriaUy this corresponds to the area of an infinitesimal square on the sphere made by moving a point on 
the surface by two orthogonal rotations (the metric is diagonal) 6 ^ 6 + d6 and 4> ^ + dcf). This method generalizes 
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to higher dimensions, in Minkowski space a sphere has the form, 

- - _ ^2 ^ 

X = iR sinip sincf) sin6 
y = iR cosip sincj) sinO 
z = iR coscf) sin9 
t = R cose 

The metric when calculated using the method just described is. 



(5.50) 



(5.51) 



^ gee 9e(t> ge^p 

9(t>e 94>4> 9m \ 
\ 9^e 9ip<t> 9il>il> 



f R^ \ 

R^sin^e 
V R^sin^e sin^cj) J 



(5.52) 



giving R^sin^O sincj) AQAcfA^j) as the surface area element. In terms of our integral over momenta we now have. 



1 



(27r)" J A 



I dip dd) dO — 7^ sin 9 sincb. 
Jo Jo Jo P 



(5.53) 

It is now a matter of grinding out the integrals by substituting mp^ completely analogous to x,y,z,t above. 

2 

The integrals with t ^ j3 are quick to do, they all vanish. For t = /? = 1, 2, 3 the integrals give — ^ and for r = /3 = 
the integral is \ so the integral is in fact equal to \g'^^ giving, 

/ (A^, (p + a) - A^, (p)) = / dS-ar (A^, {p) + a^^AA^, [p) + ...) = ^e„^^,a«fcf . (5.54) 

J (27r) (27rj J a 

The actual rest terms were, from before, 

t 
(27r) 



Rli'^ + 4^'^ = [4^ ii^f^'^ (P - ki) - ^1 - k2) - (A^. (p) -h^ k^)) (5.55) 
J (27r 



so the first term is the third term with a shift p ^ p — ki in the integrand giving —■g^eai3^ii^kfk2 and the second term 
is the same as the fourth with a shift p p — k2 giving -^eapfivk2ki, where in the latter case we have swapped 



(AX) {A2') 

fci A;2 in our formula for the shift. By asynmietry of CaiSfiv this gives, R)tv + R)iv 
consequently the anomalous/quantum axial Ward identity. 



47r' 



q^T^^x {ki,k2, q) = 2mT^^ {ki,k2,q) - -^ea/Sfi^k^k^ 

with — 4^eQ/3^iyfef ^2 being the anomaly. 

What about the rest terms in the vector Ward identity, do they cancel each other out? From before we had. 



(5.56) 



fcfVA {k,,k2,q) = if^Tr 



5 1 1 



(5.57) 
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The trace evaluates to, 

Tr [7a7^ (i*- /ki- m) 7^ 1 



m 



)] 



= Tr [7a7^ (i)- A- A2) lu {i>- A)] + Tr [to27a7^7i'] 
Use Tr [7^7.^7^] =0 

= -Tr [757A {i>- M- M)!. {i>- h)] 

Use Tr [7^ 7^,7,. 7a 7/?] = 4ie^,.a/3 



k'^T^^x{ki,k2,q) = 4 /dp - 

(27r) y \ (p - fci 



AitocfiXvk^ {p-kiY 
{p-kif 



k^r 



(p — ^2) — rn^ 



(5.58) 
(5.59) 



As before we can write this integral as a surface integral as we did in the case of the axial Ward identity. Using the 
same techniques as for the axial case we have, 



kiTfj,^x {ki,k2,q) 



(2^)" 



(p— /C2 ) — 



— 8^2 ta/Ji/Afti • 



(5.60) 



The vector Ward identity is also anomalous. 

As 1 said earlier the linear divergences essentially generate an ambiguity of the VVA amplitude, that the amplitude 
depends on how one initially defines the loop momentum. We can show this explicitly. Imagine that instead of defining 
the loop momentum as p in the initial diagram we defined it as p + a where a is some constant. For momentum 
conservation at each of the vertices this will mean that a is some linear combination of fci and k2. The amphtude for 
Tu,i,\ {ki,k2,q) is now. 



/d 
(ZTTj 

Previously the amplitude was, 

Tfj,„\ {ki,k2,q) 



P 



Tr 



y+^-»n7A7 j/_|_^_g^_m7f ji(L|_p/_^i_m7Ai 
+ 3/+<<'-m7A7 j/+9f-^-m7;i j/+^_^2-m7!' 



(5.61) 



d^p 



Tr 



1 5 1 1 

7a7 a'_</_m7i'o'_j^, _m7p 



(5.62) 



If I can simply shift the loop momentimi by the constant a, p — > p + a in the original amplitude then there is clearly 
no difference in the two amphtudes. We know that this is not the case as shifts of the variable of integration in greater 
than logarithmically divergent integrals are not trivial. Defining, 



A^i^A (p) 



Tr 



-^7A7'^^ 



y— m ''^ ' j/— 1/— m 7^^ -^—yix—ra 7/^ 

Tr[(j/+m)7A7^(j/-g'+m)7„(j/-^i+m)7f, 

(p2 — m2)((p— g)^ — rTi2)((p— fci)^— m^) 



(5.63) 



using our surface integral technology we can write. 



TfivX {ki,k2, q) - Tu,^x {ki,k2, q) = JdS'^Ur |^A^^a (p) + derivatives of A^^^a (p) ••• + 



ki 



<-^- k2 



derivatives of A^j,a (p) ■••) + 



fci k2 

[1^ V 



(5.64) 

looking at the trace the largest power of p possible in it is O (p ^), we get another three powers of p from the surface 
element p^p^. Our surface of integration is a sphere of infinite radius |p| = oo. Taking the hmit of p^ — > oo in 
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the denominator it becomes p^, consequently the integrand above is O (p^^ and the higher derivative terms above are 
O (p~ ^ ) and below. This being the case we need only concern ourselves with the first term in the Taylor expansion of 
A^^A (p + a)- Afj,^x (p), 

Ti^ux {ki,k2, q) - T^^x {h, fe, q) = - — -4 / dtp dcj) dO sin^0sin<p p^'p^ar A^^x {p) + 

(ztt) Jo Jo Jo \ 




(5.65) 



In addition, considering again the trace inside A^j^x {p) we need terms O (p^) to combine with p'^p'^ of the surface 
term and hence counter the p^ denominator so as to have the integrand non zero at \p\ = 00. We can therefore save 
ourselves several pages of trace algebra by only considering the O (p^) terms in the trace, for A/j^^x (p) this means. 



Tr[{^ + m) 7a7^ /4 + m) 7^ (^- /;i + m) 7^] Tr [inxl^ tlv 



(5.66) 



The trace is complicated and requires the use of the gamma matrix identity 7a 7/377 = Qafil-y + Q^la — ga-ylf) + 



Tr [7^ ^-f^ Jry^ ^x] 
Use identity above. 
= Pt^Tr [75 ^7^ ^7a] 

- 9^l\Tr [7^ i>j^ i> f\ 
7^7^ = 7 in last term. 

= PnTr [7^ ^7^ inx] 

- 9f^xTr [y^ ^7^ ^ ^] 
Use standard traces. 

= 4:iea„f3xP°'p'^^Pfj. 
+ ^ig^ixe-apv-iP°'p^P^ 

- ^if^oiHnXP'^PvP"' 
^nvaf}P°'p'^ type terms = 

= -^ieun-iXP^p"^ 



+ pxTr [7^ ^7i. ^7^] 

+ iea^jxP'^Tr [7^ ^j^, i>l'^l^] 

+ pxTr [7^ i>7^ i)7^] 

+ ieaij.-yxp'^Tr [in^ 

- 4.ieaf3ufiP°'p'^PX 

- '^i^apuXP'^PuP"' 



(5.67) 



Jo " dtp d(l) dO sin^esiruj) p A^^a (p) 
= /o " dV> Jo" d.^ /; dO sinHsin<p p-p^a. -^'---.gxp^- 

= '-^j^tdtpro^^I^^sinHsincP^^ 

We previously evaluated this integral in calculating the rest terms for the axial Ward identity, it was h-n'^g'^'^ 



(5.68) 



fij,vx {h, k2,q)- T^^x {h, k2,q) = 



8772 



ki ^ k2 

jJi u 



(5.69) 



Remember we said that to conserve momenta at the vertices must be a linear combination of k\ and ^2 (it is not 
also a Unear combination of g as g = fci + ^2 from simple momentum conservation), writing a = xk\ + {x — y) k2 we 
have. 



Tfj.^x (ki, k2,q) - Tfj^ux {ki, k2,q) = -^e^j^^xj {ki - fe)'' . 

does the same thing happen to the amplitude for the VVP diagram, how does it behave if the loop momentum is defined 
differently? By analogy with the VVA diagram the corresponding (surface) integrand in the case of the VVP diagram 
is the same as that of the VVA diagram but with the axial vector coupling 7a7^ replaced by the pseudoscalar coupling 

7^ 



A^. {p) = 



_ Tr [(i> + m) 75 (j)- ^ + m)j,{^- /ci+ m) 7/.] 



(p2 — m^) ^(p — q)^ — m^^ ^(p — fci)^ — m^^ 



(5.71) 
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Just as in the case of the VVA diagram the surface integral will let us take the Umit \p\ — > oo in the denominator making 
it and the integration measure will contribute p^p'^ to the numerator so Uke the VVA case we need only consider 
terms C(p^) in the trace. In other words the trace simplifies, 

Tr [(^ + m)7^(i)- j([ + m)j^{i>- M+m)j^] ^Tr[jh^ in^] = Tr [odd number of 7s] = 0. (5.72) 

Crucially the same ambiguity is therefore not true of the VVP diagram it does not change under constant shifts of the 
loop momentum p ^ p + constant, unlike the VVA diagram, the VVP amplitude is unambiguous. 

Hence had we written down our VVA diagram with p + a where previously there had been just p we would have 
had the extra term above dotted with the different 4- vectors appearing in our Ward identities! In the case of the axial 
Ward identity we would have. 



q^Tixu\{ki,k2,q) = q^Tf,^\{ki,k2,q) + ■^ef.t.x^q^ {ki - k2y 
= q^T^^x {ki, k2, q) + g^e^i/A7 (^2 ^1 ~ ^1 ^2) 



(5.73) 



Recall we calculated for the axial Ward identity. 



q^Ti^uX (fci, fc2, q) = 2mTf,^ (fci, A;2, q) - ^ea^i^^k^k^ 
q^Tnv\{ki,k2,q) = 2mTi^^{ki,k2,q) - ^ea0i^^{l + y)kfk^ 



(5.74) 



(5.75) 



For the vector Ward identity we would have, 

k^%u\{kuk2,q) = k'^Tf,,x{ki,k2,q) + ^e^,x^k'^{ki-k2y 
— k^Tfj^j^x {k\,k2,q) — ^^^ni^Xjk^k^ 

Again, in explicit calculation of the quantum vector Ward identity we found, 

fcfVA {ki,k2, q) = ^eafi,xk?kl (5.76) 



Substituting this in we find. 



k'^f^.x (fci, k2,q) = ^ea/3.A (1 - y) k-^kl (5.77) 



Hopefully it's clear that the same result is found for the other vector Ward identity, 

k^f^.x {kiM, q) = ^ea/3MA (1 - v) k^k^. (5.78) 

Therefore by a judicious choice of loop momentum in the VVA diagram (remember nothing happens to the VVP 
diagram under a constant shift of the loop momentum) we can recover the classical Ward identities for the axial vector 
and vector currents. For y = —1 we regain the classical equation of motion for the the axial current (conservation 
of the axial vector current when the loop fermions are massless) but the classical vector Ward identity is violated by 
the quantum anomaly above. For y = 1 the opposite is true and for all other values of y both the classical axial 
vector and vector current conservation equations contain extra quantum anomalies. It is a matter of convention in 
doing perturbation theory that one defines the loop momenta of such ambiguous graphs such that the vector current is 
conserved and the axial vector current contains the anomaly. 
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5.3 Fujikawa's Method for determining the Chiral Anomaly. 

In the last lecture we saw that our naive chiral Ward identity did not survive quantization. This can be attributed to the 
variance of the functional integration measure under local transformations of (only the fermion) fields giving an extra 
contribution to the path integral under the transformation which we called the anomaly function. We will now attempt 
an alternative calculation of the (Abelian) anomaly first following the method of K.Fujikawa [21-23], this is a delicate 

calculation. 

In our initial derivation of the axial vector Ward identity we found that our action of axial electrodynamics, 

S = y d^a; - ^F^^F'"' ~ ^Gi_,^G'"' +i){i^ + qy + g /A-f - m) V + Ghosts + Gauge Fixing (5.79) 

transformed under local Ua (1) transformations tp {x) — > e'^^'^^^^^^^ ip {x), 'il){x) ^ 'tp (x) e~^'^^^^^'^^ , of the fermion 
fields as, 

S + j d'^xgP {x) {di^j'X - 2imi)'y^i)) . (5.80) 

We worried that the Jacobian of the path integral measure would not remain invariant under such redefinitions of 
the fields and this is indeed the case. The essence of Fujikawa's anomaly calculation is to properly work out the 
change in the path integral measure under the transformations of the fields above. For reasons which will become 
apparent Fujikawa's anomaly calculation is only well defined in Euclidean space, we can transform the anomaly back 
to Minkowski space at the end. To first get into Euclidean space we have to make a Wick rotation to Euclidean 
space i.e. we rotate the upper index 0th component of all vectors as ia^ = and replace g^^"^ = —6^". So in 
Euclidean space we have the following identifications, ix° = x^, = 7^, do = i-^ = idi, Aq = iA^. Hence 
/D = ^^Dfjt = g^^^fjjDi^ = — 7,iDp = —j^Di — y^D2 — J^D^ — j*D4. As 7^ = 17° in Euclidean space all gamma 
matrices are anti-Hermitian (in Minkowski space only 7° was Hermitian 7^, 7^ and 7^ were anti-Hermitian). The 7^ 
matrix is still Hermitian in Euclidean space though, 

75 = ^70717^3 = 74^1^2^ ^ _^1^2^3y (5 3^) 

and we still have {7^, 7^} =0. Therefore in EucUdean space the Dirac operator above i /D = i /d-\- M + 1^ is 
Hermitian (see the end of section 1.3 for a discussion of the Hermiticity of the gauge fields). 

The next thing to do is decompose the fermion fields into eigenstates of the full Dirac operator, i.e. eigenstates of 
i jd + q y + g ~ which we will henceforth denote i p. 



i P4>i {x) = Xi(j)i {x) 
4{x)i P = \i4{x) 



+ + (5.82) 



Seeing as i /D is a Hermitian operator the eigenvalues are real! For no background fields V^, A^ these eigenstates 
are the usual Dirac eigenfunctions of definite momentum p'' as /D = /0 and A? = = w?. When the background 
fields are non-zero but well behaved and fixed, the eigenfunctions above are still the eigenfunctions ofi P but only in 
the asymptotic limit of large momentum, in which case we have p^V,A 

^ i/D<^(x)~(i^+/{/+/l7')*e-*3'- = {i> + q y +g /A^'')^e-'^-^ 

This is an important point but it is more important to note that although we will be appeaUng to this limit later we will 
at no point be assuming that we are in a plane wave basis, we are strictly dealing with the eigenfunctions of the full 
Dirac operator including the gauge field(s) i p. Note that in Euclidean space we can choose our Dirac matrices so that 
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= (p. These eigenfunctions are orthogonal and form a complete set just like in the case of the free Dirac equation, 

ih (x) = YiiaiSj (x) 

^ ^ ^ _ (5.84) 

?/; {x) = Y^icj)] {x) bi 

with fli and bj Grassmann variables ({oi,aj} = 0, {6,,6j} = 0, {ai,6j} = O), a,/? spinor indices and (/)„ (a;) 
eigenfunctions (normal numbers). The fields are now completely specified by these two (Grassmann) variables, 
consequently an integral over all possible field configurations ijj and {[> is equivalent to an integral over all possible 
combinations of values of and bi. This being the case we can now redefine the path integral measure of the fields as 

/ Di,Di, ^HfdanHf dbm. (5.86) 
So how do the Grassmann variables change under the axial transformation? 

T,jaj^j {x) — >■ (l + ip (x) 7^) Hjaj^j {x) 
^ / d'^x 4 {x) T,jaj(t)j {x) J d^x 4 (x) (l + if3 (x) 7^) ajcpj {x) (5.87) 

fli -> (Sij + J d'^x if] (x) {i\x)"/'^{x\j)) aj 
^ ai + [i J d'^x /3 {x) {i\x)'y^{x\j)) Uj. 

We wiU abbreviate the matrix J d'^x (3 {x) {i\x)'y^{x\j) by Mij so the change in ai is given by, 

ai {Sij + Mij) ttj (5.88) 

where the repeated index implies a sunnmation. Likewise we have 

i!{x)^i} {x) e+*''(^)T' 
^ ^,4>] {x) h ^ {^3^ {x) b,) (1 + (x) 7^) 

^ / d^x (e,,^] (x) &,) 0, (x) ^ / d^x (j:,4 (x) b, (1 + iP (x) 7^)) <l>i {x) (5.89) 
=> 6i ^ (<5y + / d^x iP (x) {j\x)-f^x\i)) bj 
bi^bi + i^j (/ d^xPix) {j\x)^^x\i)) bj. 

So bi transforms as, 

bi ^ {5ij + Mji) bj (5.90) 
a little different to a^. Hence the total path integral measure changes as, 

XlfdanWldb^^^jd (S, {5nj + Mnj) aj) J] / d (S, {S^j + M,„) &,) . (5.91) 

We would hke to simplify this so that the connection with the old measure is clearer i.e. we want the Jacobian whatever 
that is. To do this it will be good to know a few basic things about integrating over Grassmann variables. 

We shall now revisit our earlier discussion of Grassmann variables. Consider a function / (aj) of our Grassmann 
variable ai, we can Taylor expand it, as with functions of regular numbers, but because is Grassmann a^ai = 
—aitti = so all terms containing powers of greater than one vanish. Hence the Taylor expansion of any function 
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of a single Grassmaim variable terminates after two terms, 

/ (aO =P + a,Q (5.92) 

where P and Q are the coefficients of the expansion. 

Integration over Grassmann variables is, at least mathematically, somewhat ambiguous, consider the double deriva- 
tive of / (fli), d^.f (fli) = for any function /, so in the world of Grassmann numbers there would appear to be no 
inverse operation to left or right differentiation (for Grassmann variables differentiating from the left is different to 
differentiating from the right as da^ anticommutes with Grassmann numbers). Instead of using the formal definition 
of integration being the inverse operation of differentiation we make do with an alternative definition which makes the 
best possible analogy with that of bosonic integration. It is convention in Physics to insist that Grassmann integration 
is analogous to bosonic integration when shifting the variable of integration by a constant as this is something we often 
find ourselves doing in Physics. We want to define our integration such that under + q with Cj a constant 

number, the integral of the function over all possible values of the Grassmann variable aj is unchanged: 

Jdaif{ai) = Jd{ai + Ci) f{ai + Ci) 

= fdaif{a, + Ci) ^^^^^ 
= J da^ P + {ai + Ci)Q 
= J dui f (ai) + J dui CiQ 

To get the desired invariance under the shifts of variable above we must clearly have Grassmann integration such 
that J do, Qci = 0. To do this we define for integration of any Grassmann variable over the domain of all possible 
Grassmann numbers, 

J dtti tti = 1. 

The first of these definitions arises because we want J doj CiQ = for any function / (at) and any shift Cj. The 
second definition is merely a matter of convention, J dat is some number, it is convention to take it as one as 

a matter of normalization. These are known as the Berezin Integration Rules. Another motivation for them is that 
with these definitions the integration J da^ is equivalent to differentiating from the left (i.e. differentiation in the usual 
sense). Recall we are interested in the Jacobian involved in the transformation of a multiple integration over Grassmann 
variables. With the definitions above we have, for ci a real bosonic number. 



d{ciai) 

= (5.95) 

= / dtti tti 

d{ciai) = —dtti (5.96) 

Ci 

the Jacobian is the inverse of what it normally is! This persists to multiple integrals over Grassmann variables. Note 
that the ordering of integration of Grassmaim variables is important, defining integration of Grassmann variables the 
same as differentiation means that the integrals anticommute like the derivatives. This being the case we can get some 
useful expressions for our products of integrals, for instance we can rewrite, 

l[jdan = Ei,Ei,...Ei^-^e^i^^ -*~ j da^, j dai,... j da^^ (5.97) 

where the which are all the same up to a minus sign, which the g*i*2 - iA tensor takes care of. This gives A^! terms 
which are aU the same as JJ^ J dan taking into account the anticommuting property and the whole lot gets divided by 
A'^! giving the equivalence with Hn / '^^n- 



5.3. Fujikawa's Method for determining the Chiral Anomaly. 



77 



In our case we have, 



1 'Li . ..i^^ 



In n\ 



J dai^... / dai 



(5.98) 



where in the last Une we have taken the summation signs outside the inlegration measures they were in and used 
the rule for Grassmann integration (described above) J d(ca) = ^ J da. Because the integrals / daj^... J doj^ all 
anticommute with each other we can make the replacement, 



J daj,...J daj^ = e^-^« J da}j 



(5.99) 



hence. 



(5.100) 

Hopefully it is clear that the analysis for JJ^ J d6„ is exactly the same with the replacements a — > 6 and Mjj — > Mji 
or equivalently M — > everywhere. We can also use the fact that the determinant ofanNxN matrix is 



(5.101) 



Recall that /? (x) is infinitesimal (we are working with an infinitesimal chiral transformation) Ma <c Su so we can 
expand {da + Mii)~^ w Su — Ma. Using this definition in the working above we finally have. 



Det fe^n„/da„ 
Det(%) n„/d6„ 



(5.102) 



where 



Mi,- 



Sii - Ma i = j 



1 



Mi,. 



which is true for any matrix which can 



We can write this determinant using the identity In Det (Mi,) =lnTr 
be diagonalized by transformations of the form UMU~^. Given (3 (x) is small (we are considering an infinitesimal 
chiral transformation) so for i = j, 



]n 



Sij-i / d'^x /3{x){i\x)j^{x\j) 



i / d X P{x){i\x)Y{x\j)y i = j. 



(5.103) 



As we are interested in the trace of the above we need not concern ourselves with the case i ^ j. 

Det ^Mij) = e\p -Tr [ijd^x (3{x){i\x)j^{x\j)] 
= exp — «E„ / d'^x f3{x) {n\x}^^ {x\n} . 

Clearly we get exactly the same for fl^ J d6„. So the total transformation of the measure is, 

Y[ J dan j dhn^ (^exp - 2iE„ j d^x 13 (x) {n\x)^^{x\n)^Y\_ j ^(^n J d6„. 



(5.104) 



(5.105) 
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If we show the Dirac indices (p, a) in the exponent we have, 



^ Det {5ij -iJd'^xP {x) {j\x}j'^{x\i}) xi^ j = exp - 2iE„ / d*x (x) (1)1 ^ (x) -i%(j}nM {x) 

= exp -2lj d^X /3 {X) {j:n<Pi,p (x) <Pn,a {x)) 

Use completeness : ^,^^1^^ i^) (t>n,i3 {v) = Sa^S {x - y) 
= exp - 2i /d^a;/3(a;) (Tr [7^] .5(0)) 

(5.106) 

The sum in the integrand is badly defined, Tr [7^] = 0, (5 (0) = 00 and we need to regulate it. There are various 
methods of regularization open to us at this point, the crucial point is that they wiU aU give us the same anomaly. We 
will follow Fujikawa's regularization which was to insert a Gaussian cut-off into the sum so as to damp the contributions 
from the large eigenvalues. For a detailed discussion of all other regularization methods the reader is referred to [4]. 



'^n4>i, {x) ^^(j)n {x) = \\mM^oo^n(t>l, ip) 7^ ^Xp 

with Xn the eigenvalues of the Dirac operator i p, 

= UmM^<x)S„( n 



\2 



K ix) 



(5.107) 



7^ exp 



P 



21 



M2 



(5.108) 



This regularization is Uv (1) and IJ a (1) gauge invariant (this is easy to see if one considers the first form of the 
regulator i.e. the regulator in terms of eigenvalues instead of p). 

Now consider an eigenstate of z /D |z) with eigenvalue Aj this means that another of the eigenstates |z) = 7^ \i) 
has the opposite eigenvalue: 

i P\i) = i P^^ \i) 

= -jH P \i) . (5.109) 

= -A.|z) 

Consequently the element of the sum vanishes as all the eigenstates in the complete set with different eigenvalues are 
orthogonal. 













r x2 -1 


7^ exp 


M2 






= exp 


M2_ 



(n|n) = 



(5.110) 



The only exception occurs when the eigenvalue of the state in question is zero in which case the eigenstates need not 
be orthogonal. These so-called zero modes i P\i,0) = provide the sole contribution to the sum 



n,0 



7^ exp 



p^ 



M2 



n,0 



■ exp 



M2 



(n,0|7)^n,0= (n,0|7^|n,0> 



(5.111) 



and hence the anomaly. We can classify the zero modes according to their chirality 



7'|n+,0) 
7^|n_,0) 



Therefore our sum over all states is equal to, 

limM^(x)S„(n 



7"'' exp 



•|n+,0) 
|n_,0). 



S„ (n+,0|75|n+,0) + (n_,G|75| 
N+ - iV_ 



(5.112) 



-0) 



(5.113) 



where iV+ and are the number of positive and negative chirality zero modes respectively. The above is a statement 
of a topological theorem known as the the Atiyah-Singer index theorem [5]. The Atiyah-Singer index theorem relates 
the index of the Dirac operator to the topological charge q = — N- , this also goes under the name Pontryagin index. 
In the case of an SU (2) gauge theory it is the instanton number and for U (1) gauge theory in three dimensions it is 
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known as the (magnetic) monopole charge. For those more famiUar with topological jargon the anomaly / topological 
charge is the winding number of the map of the surface at the boundary of space to the gauge group of the theory 
in question. Consequently for trivial homotopy groups 7r„ (G) = we have no anomaly. For a detailed topological 
analysis of the anomaly the reader is encouraged to read more in [4]. 

Note that these zero modes live in a subspace of the HUbert space where effectively \i p, 7^] = 0. On the other 
hand all the non-zero modes in the rest of Hilbert space are not invariant under chiral transformations recall that 7^ |n) 
and |n) have opposite eigenvalues of P so for the space of non-zero modes \i P, 7^] = 2i Px"- It is the case that in 
our basis of eigenstates of the full Dirac operator i p, the chiral asymmetry (anomaly) is contained in the sub-space of 
zero modes. If we change to a different basis by some unitary transformation the chiral asymmetry (anomaly) will move 
into the other space of non-zero modes. We can show this for plane waves, let's Fourier transform our eigenfunctions 
so we are in a plane wave basis, 

ix) = \ -^e^'-"<^n (fc) 



limM^<x>5^n<?^l [x) 7^ exp 



= lim 
= lim 



M- 



M— >oo 



11 (fci)e-''=2-V exp 



M2 





gifci .X 







4>n (fc2)e''=i-^ 
S{h-k2). 



(5.114) 



(5.115) 



Where in the last hue we have used the completeness of 0„ , performed the sum and used the cyclic property of the 
trace in the same way as we did on the last page where we showed the Jacobian was ill-defined and needed to be 
regulated. We can rewrite the exponent using some Dirac matrix gymnastics as. 



p2 ^ ^i^^Df^D, 



(5.116) 



Where we have used the result from 1.3 for the field strength tensor F^v = [Z?^, D^]. We can also integrate over k2 
which sets fci = k2 on account of the delta function that we picked up from the completeness relation. 



'^n^i (x) 'y^(j>„ {x) = limM- 



d% 



Tr 



e *'=i-'^7^ exp 



M2 



M2 



(5.117) 



The next thing to do is move the factor of e**^! through to meet the other one which gives a delta function of ki and 
k2. This is easy but we need to look out for the differential operator in D^D^ acting on the e**^! "^ factor. In the limit 
M — > 00 we can 



e »'=i-^75 exp 



M2 



Ak\.x 



5„-iki.x ( D'^Df^ 



1 /Dt'D, 



■ 7 6 



M2 



M2 



Aki.x 



(5.118) 



It is now necessary to simplify the D^D^^ term, all the while we have to consider that this is an operator equation, i.e. 
that D^D'^ is acting on something else on the left besides 6**^^ 



ki.x (^fl j_ A "Sf^fllJ-J- AIJ'\^iki-x 



(5.119) 



Now we move e''^^-^ left through the next Z)^, this has the same effect as above and in addition the two exponentials 
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cancel, 



g-ifcl.X£)^£)/.gi/ci.X 



{D^ - ki^) {D^ - k1) . 



(5.120) 



Repeating this for the following terms in the series e '^i-^i ( ^^f?" ) e''^^'^ gives a completely analogous result i.e. 
simply replace D ^ D — k\ and remove the exponentials. This being the case we have, 



7''' exp 



= lim 



M^oo 



"Tr 



(D''-fci')(D^-fci^) ^[7^7nJ'M 



7^ exp 



limM^oo / A^^M^e+'^i'^^-Tr 



M2 



+ 



M2 



7^ exp 



M 



+ 



M2 



+ 



M2 



2M2 



(5.121) 



where in the last step we simply rescaled the momentum integral as ki — > Mki and used the fact that F^^^, is n u 
symmetric. 

In Taylor expanding the exponential the first term that will be non-zero is the one quadratic in F^i,. This is 
because of the 7^ in the trace (the trace was over the eigenfunctions of /D in position space (i.e. x) but also naturally 
over the Dirac indices). The first non vanishing trace including a 7^ is the trace of 7^ with four other 7 matrices, 
Tr [75 

^A'^i'^a^/^j — _4gA'i'a/3 (jjj Euclidcan space). The first term in the expansion to have this property which is 
least suppressed by powers of M is clearly, 

1 (l^'YF^.^ ' 



2M2 



(5.122) 



This is also the only non vanishing term in the M ^ 00, it is the only such term with four powers of M in the numerator 
and the denominator (we picked up a factor M'^ on rescaling the momentum), all other terms with non- vanishing traces 
will have at least one power of M in the denominator and so vanish as M ^ 00. Therefore, 



Y.ncj>i {x) j^<Pn {x) = limM^oo / ^M^e+^'^^-^Tr 7^ ^ (2^^^) ' 
= / g^e+'^i '^^MTr [75 i {l^YF^^f 



(5.123) 



(the minus sign appears from the trace over Dirac indices in Euclidean space Tr [7^7^7''7"7^] = — 4gA'i'a/5) Now we 
recall that in going from Minkowski space to Euclidean space we have /cf fci^ = —ki^ki^ where the fci^/ci^ vectors 
are Euclidean. We can therefore perform the integral over ki_ {F^^ does not depend on k\) as it is now just a product 
for four regular Gaussian integrals (one for each dimension). 



/ 



d^fci e-'^i"^!" = TT^ 



E„(n |7^| n) = -J-^e^^"^F^,F^p 



327r2 

Going back a few pages we found that the path integral measure changed as. 



(5.124) 
(5.125) 



H J dan J dbn ^ (^exp - 2iE„ J d^x p {x) {n\x)^^{x\n)^ 11 /'^"" / (^-126) 



substituting in the result above we have. 



n„/da„/d6„ ^ {^xpji^Jd'xP{x)e^^-^-^F^^F„p)UnIdan!dbn 
= (exp -Jd^x(3 {x) A [A^]) n„ / da„ / d6„ 



(5.127) 
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where, 

A [A^] = -^e^^'^^^F^^F^p (5.128) 

is the so-called singlet / axial / chiral anomaly. To get the anomaly in Minkowski space we need to make one change 
which is that taking the trace of the 7 matrices above in Minkowski space gives —Aie^""^ Fni,Fai3, which means it is 
different by a factor —i, 

A [A^\ = ^e'^^'^^F^^F^p. (5.129) 

So in performing the usual transformation of variables technique to derive the axial current conservation equation we 
have the same terms as before but now we also have the anomaly term, 

^ {d^fj,) = 2im(V;7'V') + {A [A^]). (5.130) 

It turns out that this expression is the same as that which we calculated in perturbation theory from the triangle diagrams, 
to do this you consider the field strength tensor F^^, associated with the two photons in the triangle diagram with 
momenta k\, and polarization vectors ei and €2. This is one of a long list of peculiarities of the anomaly, that the 
lowest order calculation of the anomaly {i.e. with the triangle diagrams) receives no higher order corrections, it agrees 
with the non-perturbative calculation in this section. 

To finish off we list briefly some of the broader issues and spin offs associated with this calculation, unfortunately 
time and space does not permit more detail. 

The first thing is related to an old problem associated with the electromagnetic decay of the pion tt" 77, this 
decay is a triangle diagram of the form calculated in the previous section. Pions are pseudoscalar particles and in the 
limit of soft interactions {q « 0) we can effectively consider them as a point like particle with a pseudoscalar coupling 
{i.e. only a 7^ at the vertex) thus their decay basicaUy looks like the VVP diagram in section 5.2. However the naive 
Ward identity we derived between the VVA and VVP diagrams told us that, 

q^T^i^\{ki,k2,q) =2mTfj,,y{ki,k2,q) (5.131) 

i.e. that the amplitude of the VVP diagram is effectively zero in the limit g ^ 0, or equivalently, the 7r° doesn't decay! 
This was called the Sutherland-Veltman Paradox. It is one of the great successes of the discovery of the anomaly that 
the 7r° lifetime agrees very well with experiment when one takes into account the anomaly in the triangle diagrams of 
low energy effective theory of the tt'^. Such a calculation of the 7r° lifetime can be found in [3]. 

Anomalies are more well known these days for their role in constraining physical theories. The constraints arise 
because innocuous as it may seem the fact that anomalies violate our classical Ward identities is another way of 
saying that the gauge invariance associated with the anomaly does not exist at the quantum level. The existence of the 
quantum Ward identities is crucial to the renormalizability of a theory, we saw how they killed off the nasty quadratic 
divergences in QED for example. It also enabled us to derive relations between divergent diagrams which in turn 
gave a cancellation of those divergences (perhaps a better way of saying this is that it gave us relations between the 
renormalization coefficients). The fact that these identities are not true ruins the Ward identities and so we lose the 
abiUty to have such cancellations, this was very clear in the case of trying to do perturbation theory with the VVA 
diagram which showed an ambiguous amphtude, the naive Ward identity between the VVA and VVP diagrams was 
untrue. It is also the case that anomalies give more nonsense in so far as the S-Matrix acquires a gauge dependence, 
see [9]. A sensible unambiguous (and renormalizable) theory is necessarily anomaly free, this is the aforementioned 
constraint which aids model building. Had we considered a non-Abelian gauge theory in the treatment of the last two 
sections we would have had various factors involving the gauge group generators appearing, around the vertices of the 
triangle diagrams in 5.2 and in the transformations of the fermion fields in 5.3. The modification to the calculations 
is trivial and the calculations proceed as they do in 5.2 and 5.3 (but with more mess). The upshot is that the anomaly 
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contains a factor proportional to the gauge group generators, 

where the trace is over the indices of the internal symmetry space. Substituting in the generators of Electroweak theory 
the anomaly from the VVA diagrams can yield a term such that, 

tr [e [th'' + tH'']] = tr [e + Yt'']] ex 5''hr Q. (5.133) 

Thus anomaly cancellation in the Standard Model requires that the sum of the electric charges of all the fermions 
should vanish^. 



^Remember a colour factor Nc = 3 when adding up the quarks charges. 
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